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1.0 Introduction 


The practice of limiting the shaker force in vibration tests was instigated at the NASA Jet 
Propulsion Laboratory (JPL) in 1990 after the mechanical failure of an aerospace component 
during a vibration test. Now force limiting is used in almost every major vibration test at JPL 
and in many vibration tests at NASA Goddard Space Flight Center (GSFC) and at many 
aerospace contractors. The basic ideas behind force limiting have been in the literature for 
several decades, but the piezo-electric force transducers necessary to conveniently implement 
force limiting have been available only in the last decade. In 1993, funding was obtained from 
the NASA headquarters Office of Chief Engineer to develop and document the technology 
needed to establish force limited vibration testing as a standard approach available to all NASA 
centers and aerospace contractors. This monograph is the final report on that effort and 
discusses the history, theory, and applications of the method in some detail. To facilitate the 
application of the method, a more concise description of the key aspects of the approach is 
presented in a complementary guidelines document, which is available separately and is also 
contained herein as Appendix C. 
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2.0 History 

Many pioneers of aerospace mechanics, as early as the 1950’s, recognized that ignoring 
the low mechanical impedance of lightweight mounting 

unrealistically severe vibration tests of aerospace equipment Unfortunately, J hei ^ 
wakings were not heeded until recently, primarily because the tools wen: not available 
to conveniently simulate mechanical impedance in vibration tests. Also, more realistic 
“SS essential in today's faster-better-cheaper environment, whtch is 
incompatible with the traditional approaches of high design margins and extensive 
developmental testing. 

Some key background references and their contributions to the force limited vibration 
technology described in this monograph are summarized in chronological order in this 
sechon ^ These and other supportini references are included in the bibliography of 

Appendix B. 

nt k. |qc<; r i n Hesrribes the problem of overtesting at resonances of the test item 
whth results om theSdXact.ce of enveloping the peaks in the field acceleration 
" iiafda a for b “h vibration and shock tests. He proposed a complex, conceptual 
E to which the impedance of the mounting structure would be simultaneously 
measured with a small shaker and emulated by the test shaker. 

Morrow 1960 [2] warns against ignoring mounting structure impedance m both 

vibration and shock tests and points out that impedance concepts fjT 1 ™ “ Xdie 
engineers are largely unknown to mechanical engineers. He describes exact impeaan 
rSaSsin?force ttansducers between the shaker and test djm 
difficulties in simulating impedance exactly, phase included. (It is still impractical to 
specify and control the shaker impedance exactly.) 

Salter 1964 131 calls for two test improvements to alleviate overtesting. 1) multi-point 
control toreduce the impact of fixture resonances and 2) force limiting to account for 
Ok vibrMion absorber effect at test item resonances. He proposes a very simple method 
Of force limit, i.e. the force is limited to 1 .5 times the mass times the 

peakTccderation, i.e. the accelerate, specification. His 

~ review of the force data obtained in the system acoustic tests of the Cassini spacec 
provides the impetus for what in this monograph is called the semi-empirical method 

predicting force limits. 

Ratz 1966 Ml who was the chief engineer for MB electronics, designs and tests a 
new shaker equXer which uses force feedback to simulate the mechanical impedance 
of die equipment mounting structure (foundation). The last sentence of his paper, Use 
of the new equalizer, therefore, can make a dramatic improvement in vibration 
simulation and would seem to be the harbinger of a significant advance in the state o 
the art of vibration testing,” looks forward to what is just now becoming a reali y. 

Heinricks 1967 [5] and McCaa and Matrullo 1967 [6] in complementary papers 
describe an analysis and test, respectively, of a lifting body re-entry vehicle using fore 
limiting to notcha random vibration acceleration spectrum. A complete modal model 
ndudmg tte effilta mass concepts discussed herein, are developed in the analysis^ 
The analysis also includes a comprehensive finite element model (FEM) simulation o 
SSSl. test, in which the test input forces arc limited to the sttuctural 
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mit-load criteria. A vibration test of a scale model vehicle is conducted using single- 
ax s impedance-head force transducers to measure the total force input, and the § 
notching is implemented manually based on force data from low level tests. 

Painter 1 967 [7] conducts an experimental investigation of the sinusoidal vibration 
intirf? of f aircraft components using both force and acceleration specifications The 
nterface forces and accelerations between simulated equipment and an aircraft fuselage 
were measured and enveloped, and these envelopes are used to control shaker g 
sinusoidal vibration tests of the equipment. It is found that the procedure largely 
eliminated the high levels of overtesting introduced by the conventional approach. 

® 19 f 7 ^ ^ nd Rubin 1967 PI in two complementary papers translate and 

extend the electrical engineering impedance concepts into mechanical engineering terms 
O Hara points out an important distinction between impedance, the 
applied velocity and mobility, the ratio of velocity to an applied fome ?n the 

lort^/h 'h Cnt ° f , jy ed ® nce » for multiple drive points or degrees of freedom, the motion 
r each degree of freedom other than the one of interest must be blocked i e 
constramed to have zero motion, which is difficult to realize experimentally ' The 
mobility concept, on the other hand, requires simply that the force at other degrees of 
freedom be zero i.e. free boundary conditions. The impedance and moSiSSStoto 

m Ca / S ' In the , sim P le . case where th e impedance matrice is diagonal (i e the 
° f the force to the motion for different degrees of freedom are zero), this 
mction does not exist, and the impedance matrix components are simply the 
reciprocal of the corresponding mobility matrix component and equal to whaf O’Hara 
ca^ls pseudo-impedance. (In this monograph, an impedance-like quantity the apparent 

rnmniw he ratl o °f force to applied acceleration, is most often employed and to avoid 
complications, which are considerable, the off-diagonal terms of the impedance matrix 

usefnl P tor y ne S lected J Rubln developed transmission-matrix concepts P which are very 
useful for coupling systems together and for analyzing vibration isolation. ^ 

Murfin 1968 [10] develops the concept of dual control, the first of several at Sandia 

,0 the lech k n0 ' o » r of force limiting, which i., one of 
tne Key elements ot the force limiting approach described in this monograph He 

proposes that a force specification be developed and applied in a manned completely 
analogous to the acceleration specification. The force specification is to be thesmoothed 
envelope of the force peaks in the coupled system res^nse, i.e. the field! enviro^nT 
Recognizing that such field force data are not readily available, he proposes a melhod of 
enving the forces from the product of the acceleration specification and the smoothed 
apparent mass of the test item. He ignores the mounting structure impedance This 
approach is very similar to that used in the semi -empirical method described in this 
"!°" 0 r gra P l ?' In the dual con j ro1 concept, the vibration controller does not give 
preferential treatment to either the acceleration or the force specification. The force and 
cceleration control signals are analyzed in narrow bands, and in each band the shaker 

ivniruM^ 11 ^ ls . ad J usted until one or the other of the specifications is attained. This 
yp y results in acceleration control off resonances and force control and 
acceleration notches, at the resonances. Many older vibration test controllers have an 
extremal or peak control mode in which the largest of several control signals are 
compared to the acceleration specification, but do not allow separate reference 
speafications for acceleration and force or response. To use these older controllers for 

rink ’l l nccessary t( ? use a shaping filter to condition the force or response 
contro 1 s^nul Most newer vibration test controllers have the capability to implement 

“ P „r syss w,,h some minor variati °" s 
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Scharton 1969 [11] develops special, multi-modal, vibration test fixtures which had 
enhanced modal densities and low rigidities, to mechanically simulate the impedance of 
large flight mounting structures. 

Witte 1970 [12] proposes a method of controlling the product of the force and 
acceleration which is applicable when no information is available on either the test item 
or the mounting structure mechanical impedance. 

Witte and Rodman 1970 [13] and Hunter and Otts 1972 [14] continue to 
pursue the calculation of the force specification by multiplying the acceleration 
specification times the smoothed apparent mass of the test item; like Murfin and their 
colleagues at Sandia, they ignore the mounting structure impedance. They use simple 
parametric models to interpret field data and to study the dynamic absorber effect of the 
payload at resonance, and they develop special methods of smoothing the test item 
apparent mass. 

Wada, Bamford, and Garba 1972 [15] develop a technique for obtaining an 
equivalent single-degree-of-freedom system (SDFS) for each eigen-vector when the 
dynamic characteristics of the structure are available in the form of a finite element 
model (FEM) or as test data. The masses of these SDFS’s are called the “effective 
mass” which may be defined as the mass terms in a modal expansion of the drive point 
apparent mass of a kinematically supported system. The effective masses are 
independent of the modal normalization, and the sum of the effective masses of all 
modes is equal to the total mass. A complementary term is the residual mass, which is 
defined as the difference between the total mass and the effective masses of all the 
modes with natural frequencies lower than the frequency of interest. It is a consequence 
of its definition that the residual mass must be a decreasing function of frequency 
(Foster’s Theorem). The effective mass is very important concept, because it provides a 
way to quantify, from either FEM’s or test data, the mass of a structure as a function ot 
frequency. The effective mass and residual mass are used in this monograph to 
characterize the mechanical impedance of both the mounting structure (the source) and 
the test item (the load). 

Martini 1983 [16] describes the advent of the piezo-electric, quartz, multi- 
component force transducer, which is certainly the most important enabling factor in 
making force limited vibration testing a reality. The development of quartz multi- 
component force transducers started in 1965 with a Swiss government project to 
provide a very stiff sensor to measure cutting forces on machine tools, continued in the 
seventies with the development of biomedical applications, and finally came into its 
own in the early 80’ s when the automotive industry began using six component 
dynamometers for measuring tire and wheel loads. Piezo-electric, quartz transducers 
offer the following crucial advantages over strain transducers, which have traditionally 
been used to measure force: 1 . extreme rigidity and therefore high resonance frequency, 
2. wide dynamic range ~10 4 , 3. large span-to-threshold value ~10 , 4. natural 
resolution into orthogonal components, 4. compactness, 5. wide temperature range - 
200C to +200C, 6. low cross-talk, and 7. direct measurement of total force, rather than 
strain. As seen from the previous references, many of the ideas behind mechanical 
impedance simulation, and more specifically force limiting, were reported in the 
literature some thirty years ago. Yet force measurements were seldom made in vibration 
tests until the 1990’s, when quartz multi-component force transducers became readily 
available, primarily as a result of their development for other markets. 

Judkins and Ranaudo 1987 [17] conduct a definitive series of tests quantifying the 
degree of overtesting in conventional aerospace vibration tests. Their objective is to 
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compare the damage potential of an acoustic test and a conventional random vibration 
test on a shaker. The test item consists of a three slice mock-up of a RF component 
weighing approximately 15 lb. and containing three simulated circuit boards. In the 
acoustic test, the component is mounted on a 0.5” thick honeycomb panel with 
dimensions of 28.85” by 45.22”.The study shows that the shaker resulted in an 
overtest factor (ratio of shaker to acoustic test results) of 10 to 100 for peak spectral 
densities and a factor of ten for G rms’s. They point out that significant savings in 
design schedules and component costs will result from reduced vibration test levels 
which are developed by taking into account the compliance of the mounting structure in 
the vibration tests of spacecraft components. 

Sweitzer 1987 [18] develop a very simple method of correcting for mechanical 
impedance effects during vibration tests of typical avionics electronic equipment. In 
essence, the method is to let the test item have a resonance amplification factor of only 
the square root of Q, rather than Q as it would on a rigid foundation. This is 
implemented in the test by notching the input acceleration by the same factor, i.e. the 
square root of Q. This simple method is attractive because it requires no additional 
instrumentation or knowledge of the actual impedance of either the test item or the 
mounting structure. 

Piersol, White, Wilby, Hipoj, and Wilby 1988 [19] conduct a definitive study 
of the causes and remedies for vibration overtesting in conjunction with Space Shuttle 
Sidewall mounted components. (It was while working as a consultant on this program 
that the author became interested in force limited vibration testing.) They compile an 
extensive biography, contained in Appendix A of their Phase I Report, of literature 
relating to impedance simulation and the vibration overtesting problem. One aspect of 
their study is to obtain impedance measurements on the shuttle sidewall and correlate 
the data with FEM and semi-empirical models. They propose as a force limit the 
blocked force , which is the force that the field mounting structure and excitation 
would deliver to a rigid, infinite impedance, load. Unfortunately, results show that the 
blocked force is still very conservative for most aerospace applications, so that the 
overtesting problem is not much alleviated with this particular choice of a limit. A 
further drawback of the blocked force is that it does not take into account the impedance 
of the test item, which is readily available in a vibration test incorporating force 
transducers. 

Scharton and Kern 1988 [20] propose a dual control vibration test in which both 
the interface acceleration and force are measured and controlled. They derive an exact 
dual control equation which relates the interface acceleration and force to the free 
acceleration and blocked force. The exact relationship is of little practical value because 
current vibration test controllers cannot deal with phase, and the exact characteristics of 
the mounting structure impedance can not be determined. Alternately, they propose an 
approximate relationship, for dual extremal control, in which the exact relation is 
replaced by extremal control of the interface acceleration to its specification and the 
interface force to its specification, as discussed by Murfin. 

Scharton, Boatman, and Kern 1989 [21] describe a dual controlled vibration test 
of aerospace hardware, a camera for the Mars Observer spacecraft, using piezo-electric 
force transducers to measure and notch the input acceleration in real time. Since the 
controller would not allow a separate specification for limiting the force, it was 
necessary to use a shaping filter to convert the force signal into a pseudo-acceleration. 
One of the lessons learned from this project was that the weight of the fixture above the 
force transducers should represent a small fraction (less than 10%) of the test item 
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weight. Otherwise above approximately the first mode of the test item, the force signal 
will be dominated by the force required to vibrate the dead weight of the fixture. 

Smallwood 1989 [22] conducts an analytical study of a vibration test method using 
extremal control of acceleration and force. He finds that the method limited the 
acceleration input at frequencies where the test item responses tend to be unrealistically 
large, but that the application of the method is not straightforward and requires some 
care. He concluded that the revival of test methods using force is appropriate 
considering the advances in testing technology in the last fifteen years, and that the 
method reviewed shows real merit and should be investigated further. 

Scharton 1990 [23] analyzes dual control of vibration testing using a simple two- 
degree-of-freedom system. The study indicates that dual controlled vibration testing 
alleviates overtesting, but that the blocked force is not always appropriate for the force 
specification. An alternative method is developed for predicting a force limit, based on 
random vibration parametric results for a coupled oscillator system described in the 
literature. 

Smallwood 1990 [24] establishes a procedure to derive an extremal control vibration 
test based on acceleration and force which can be applied to a wide variety of test items. 
This procedure provides a specific, justifiable way to notch the input based on a force 

limit. 

Scharton 1993 [25] describes application of force limited vibration testing to nine 
JPL flight hardware projects, one of which is the complete TOPEX spacecraft tested at 
NASA Goddard Space Flight Center. Two of the cases include validation data which 
show that the force limited vibration test of the components are still conservative 
compared with the input data obtained from vibration tests and acoustic tests at higher 

levels of assembly. 

Scharton 1994 [26] describes two applications of force limiting: the first to the Wide 
Field Planetary Camera II for the first Hubble telescope servicing mission, and the 
second to an instrument on the Cassini spacecraft. 

Scharton 1995 [27] devises a method of calculating force limits by evaluating the test 
item dynamic mass at the coupled system resonance frequencies. Application of the 
method to a simple and to a complex coupled oscillator system yields non-dimensional 
analytical results which may be used to calculate limits for future force limited vibration 
tests. The analysis for the simple system provides an exact, closed form result tor the 
peak force of the coupled system and for the notch depth in the vibration test. For 
example, using the simple system results with Q=50 and equal impedance of the flight 
mounting structure and test item, the input acceleration will be notched by a factor of 
31.25 relative to a conventional test. The analysis for the complex system provides 
parametric results which contain both the effective modal and residual masses of the test 
item and mounting structure and is therefore well suited for use with FEM models. 

Scharton and Chang 1997 [28] describe the force limited vibration test of the 
Cassini spacecraft conducted in November of 1996. Over a hundred acceleration 
responses were monitored in the spacecraft vibration test, but only the total axial force 
is used in the control loop to notch the input acceleration. The force limit specified in 
the spacecraft vibration test plan is used in the test without any modifications, and many 
of the major equipment items on the spacecraft reached their flight limit load. The force 
limit for the complete spacecraft vibration test, as well as the limits for many of the 
Cassini instrument vibration tests, are developed using a simple, semi-empirical method 


2-5 


which requires only the acceleration specification and data from a low level pre-test to 
determine the apparent mass of the test item. This semi-empirical method of predicting 
force limits is validated for the instrument tests by comparisons with two-degree-of- 
freedom analytical models and with interface force data measured at the 
instrument/spacecraft interface in acoustic tests of the Cassini spacecraft DTM structure 
The instrument force limits derived with the semi-empirical method are generally equal 
to or less than those derived with the two-degree-of-freedom method, but are still 
conservative with respect to the interface force data measured in the acoustic test. 
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3.0 Structural Impedance 


3.1 The Vibration Overtesting Problem 

There are historically three solutions to the vibration overtesting problem. 1) build it like a 
brick”, 2) mechanical impedance simulation, and 3) response limiting. 

Some aerospace components are still “built like a brick” and therefore can survive vibration 
overtesting and perhaps even an iterative test failure, rework, and retesting scenario. In a few 
cLes tlds may even be the cheapest way to go, but the frequency of such cases is certain y 
much less than it used to be. The two historical methods of alleviating overtesting, impedance 
simulation and response limiting, are both closely related to force limiting. 

3.1.1 Impedance Simulation 

In the 1960’s personnel at NASA’s Marshall Space Flight Center (MSFC) developed a 
mechanical impedance simulation technique called the “N plus one stnicture” concept, which 
involved incorporating a portion of the mounting structure into the vibration test A common 
example would be the vibration test of an electronic board, mounted in a black tox In 
addition acoustic tests were often conducted with the test items attached to a flight 
mounting structure. In all these approaches where a portion of the mounting s ^ture, °r 
simulated mounting structure, is used as the vibration test fixture, it is preferred that the 
acceleration input be specified and monitored internally at 

structure and the test item. If instead, the acceleration is specified externally at the interface 
between the shaker and mounting fixture, the impedance simulation benefit is greatly 
decreased. (When the input is defined internally, the “N plus one structure approach is similar 
to the response limiting approach, discussed in the next section.) 

A second example of mechanical impedance simulation is the multi-modal vibration test fixture 
rill which was designed to have many vibration modes to emulate a large flight mounting 
structure. This novel approach was used in one government program, a Mariner spacecraft but 
fell bv the wayside along with other mechanical impedance simulation approaches as being too 
Sw and expensive. In addrtion, the concept went against the convent, onal wtsdom 
of making fixtures as rigid as possible to avoid resonances. 

Mechanical impedance simulation approaches are seldom employed because they require 
additional hardware and therefore added expense. Two exceptions which may find acceptance 
in this new low cost environment are: 1) deferring component testing until t 
assembly e.g. the system test, when more of the mounting structure is automatically present, 
and 2) replacement of equipment random vibration tests with acoustic tests of the equipment 
mounted on a flight-like plate, e.g. honeycomb. 


3.1.2 Response Limiting 

Most institutions have in the past resorted to some form of response limiting as a means of 
alleviating vibration overtesting. Response limiting is analogous to force hunting, but gener y 
more compifcS™ d dependent on Lalys, s. Response limiting was used for several decades 
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at JPL but has now been largely replaced by force limiting. In both response and force limiting 
the approach is to predict the in-flight response (force) at one or more critical locations on the 
item to be tested, and then to measure that response (force) in the vibration test and to reduce 
or notch, the acceleration input in the test at particular frequencies, so as to keep the measured 
response (force) equal to or below that limit. F 

In the case of response limiting, the in-flight response is usually predicted with FEMs. This 
means that the prediction of response requires an FEM model of the test item and the directly 
excited supporting structure and the same model is typically used to design and analyze the 3 
oads in the test item. In this case the role of the test as an independent verification of the design 
and analysis is severely compromised. In addition, the model has to be very detailed in order to 
predict the m-flight response at critical locations, so the accuracy of the predictions is usually 
suspect, particularly at the higher frequencies of random vibration. By contrast the interface 
lorce between the support structure and test item can be predicted with more confidence and 
depends less or not at all on the FEM of the test item. 

In addition, it is often complicated or impossible to measure the responses at critical locations 
on the vibration test item. Sometimes the critical locations are not accessible, as in the case of 
optical and cold components. In the case of large test items, there may be many response 
locations of interest; hundreds of response locations may be measured in a typical spacecraft 
test_ For this reason, some institutions rely completely on analysis to predict the responses in 

m i he tCSt ’ and then a pnori shape the input acceleration for the test in order to equate 
the flight and test responses. Since the uncertainty in the predictions of the resonance 
frequencies of the item on the shaker is typically 10 to 20 %, any notches based on pre-test 

analysis must be very wide, and may result in undertesting at frequencies other than at 
resonances. 

There is one form of response limiting which is conceptually identical to force limiting i.e. 
uniting the acceleration of the center-of-gravity (CG) or mass centroid of the test item. By 
Newton s second law, the acceleration of the CG is equal to the external force applied to the 
body divided by the total mass. It is indeed much easier to predict the in flight responses of the 
test item CG, than the responses at other locations. The CG response is typically predicted with 
F n ^ S u Slng ° n y 3 umped mass to represent the test item. At JPL a semi-empirical curve 
called the mass-acceleration curve is usually developed early in the design process to predict the 
CG response of payloads. Also, any method used to predict the in-flight interface force 
obviously predicts the CG acceleration as well. 


The problem with CG response limiting in the past has been that it is difficult or impossible to 
measure the acceleration of the CG with accelerometers in a vibration test. Sometimes the CG 
is inaccessible, or there is no physical structure at the CG location on which to mount an 
accelerometer. However, there is a more serious problem. The CG is only fixed relative to the 
structure, when the structure is a rigid body. Once resonances and deformations occur it is 
impossible to measure the CG acceleration with an accelerometer. Furthermore attempts to 
measure the CG response usually overestimate the CG response at resonances, so limiting 
based on these measurements will result in an undertest. However, the CG acceleration is 
uniquely determined by dividing an interface force measurement by the total mass of the test 
item, this technique is very useful and will be discussed subsequently in conjunction with 
quasi-static design load verification. 


3.1.3 Enveloping Tradition 

The primary cause of vibration overtesting is associated with the traditional, and necessary 
practice of enveloping acceleration spectra to generate a vibration test specification. In the past 
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3.1.4 Dynamic Absorber Effect 

The dynamic absorber effect [30] may be explained 
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Acceleration Input to Electronic Boxes (g2/Hz) 



FIGURE 1. Measurements of Random Vibration Acceleration Spectra on Honeycomb Panel 
Near Mounting Feet of Electronic Boxes in TOPEX Spacecraft Acoustic Test 
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FIGURE 2. Simple Two-Degree-of-Freedom System (TDFS)directly excited oscillator When 
he damping of the secondary oscillator is small rather than zero, the motion of 
the primary mass is small rather than zero. 


3-4 









force. Fig. 3b is the magnitude of the coupled system interface acceleration, and Fig. 3c is the 
magnitude of the load apparent mass. 

In Fig. 3, notice first that the interface force and interface acceleration both have peaks at the 
two coupled system resonance frequencies of 0.62 f 0 and 1.62 f 0 . Notice further that the 
interface acceleration has a notch of depth Q 1 at the load fixed-base resonance frequency f 0 , 
where and the load apparent mass has a peak of height Q. This notch in the interface 
acceleration is just the dynamic absorber effect. This example illustrates the general dynami 
absorber result, that the frequency spectrum of acceleration at the interface between 
spacecraft and instrument will have notches at the fixed-base resonance frequencies ot the 

instrument. 

The example in Fig. 3 may also be used to illustrate the overtesting resulting from enveloping 
the interface acceleration/and how force limiting will alleviate this overtesting. In ^ coupled^ 
system, the interface force peak of 80 at the lower resonance frequency of 0.62 f results trom 
multiplying the interface acceleration peak of 50 by the load apparent mass value of 1.6. 
conventional vibration test without force limiting, the corresponding shaker force wouldb o 
2500 which is the interface acceleration envelope of 50 times the load apparent mass p 
50 at the load resonance frequency f. With force limiting, the input acceleration wouldbe 
notched at the load resonance frequency f 0 to reduce the shaker force by a factor of 2500/8 

31.25. 


3.2 Dual Control of Acceleration and Force 

Conventional vibration tests are conducted by controlling only the acceleration input to the test 
item In theory, if the frequency spectrum of the acceleration input in the test, including peaks 
and valleys, were identical to that of the interface acceleration in the flight mounting 
configuration, and if the boundary conditions for other degrees-of-freedom (rotations, etc.) 
were the same as in the flight configuration, then the interface forces and all the responses 
would be the same in the test as in flight. However, this is seldom the ^ se ’ USe 
of the necessity of using a smoothed or an enveloped representation of the flight interface 
acceleration as the test input, and secondarily because of frequency shifts ^s^i^ed with e 
unrealistic restraint of other degrees-of-freedom by the shaker mounting. It has been found tha 
the dual control of the acceleration and force input from the shaker al eviates the overtesting 
problem associated with conventional vibration tests using only acceleration contro . 

3.2.1 Thevinen and Norton’s Equivalent Circuit Theorems 

Consider a source, consisting of a voltage source in series with a source impedance, which is 
connected to a load [31]. If we adopt the mechanical analogy m which force is current and 
velocity is voltage (Unfortunately, this is the dual of the mechanical analogy used in 
reference^ 1].), Thevinen’ s equivalent circuit theorem may be stated in mechanical terms as. 


A = A 0 - F/ M, 


( 1 ) 


where A is the source-load interface acceleration, A 0 is the free acceleration (i.e., th 
acceleration that would exist at the interface if the load were removed), F is the interface force, 
and M is the source apparent mass measured at the interface. (Apparent mass is discussed i 
the next section.) All the terms in Equation 1 are complex and a function of frequency. 
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Using the same electo-mechanical analogy, consider a source, consisting of a current source in 
parallel with a source impedance, connected to a load. Norton’s equivalent circuit theorem 
stated in mechanical terms is : 


F = F -AM, 


( 2 ) 


where F is the blocked force (i.e., the force that would be required at the interface to make the 


motion zero). 


Equations (Eq.)l and 2 may be manipulated to eliminate both F and A yielding the following 
relationship between the blocked force, free acceleration, and the source apparent mass. 


F / A = M- 


(3) 


3.2.2 Dual Control Equations 

Alternately, the source apparent mass M may be eliminated from Eqs. 1 and 2, to yield the 
following [20]: 

1= A/A 0 +F/F 0 , ^ 

which provides a theoretical basis for dual control of vibration tests. 

Equation 4 is exact but difficult to apply because the terms on the right hand side are complex 
and complicated functions of frequency. The phase of the inputs and the impedance are dif 
to determine analytically or experimentally, although some exploratory work on this problem 
was conducted some 25 years ago [4,14], Little recent work on exact mechanical impedance 
simulation is available, and most commercially available vibration test controllers can not 
control phase angle to a specification. 

An alternative, approximate formulation for the control of vibration tests is provided by the 
following extremal Eqs. [7]: 


IAI/IAJ < 1 and IFI/IFJ < 1, 


(5) 


in which A represents the acceleration specification and F s represents the force specification, n 
Eq. 5 the free acceleration and blocked force of Eq. 4 are replaced by the corresponding 
specifications which envelope the interface acceleration and force in the coupled system. With 
extremal control, the shaker current is adjusted in each narrow frequency band so that the larger 
of the two ratios in Eq. 5 is equal to unity. At frequencies other than the test item resonances, 
the acceleration specification usually controls the test level; at the resonances, the base reaction 
force increases and the force specification limits the input. 

Most vibration controllers have the capability for extremal control, but older controllers allow 
only one reference specification. To implement dual control in this case, a filter must be used to 
scale the shaker force feedback signal to an equivalent acceleration [17].New controllers allow 
separate specifications for limit channels, so Eq. 2 may be directly implemented. Force limiting 
has been used primarily for random vibration tests, but the application to swept sine tests is 
also practical and beneficial. 


3-7 


3.3 Structural Impedance Characterization 


3.3.1 Apparent Mass 


In this monograph structural impedance will be characterized as “apparent mass”, which is the 
preferred name for the frequency response function (FRF) consisting of the ratio of reaction 
force to prescribed acceleration [32]. (Apparent mass symbols will be underscored in this 
monograph to distinguish them from other mass quantities.) The force and prescribed 
acceleration in the apparent mass usually refer to the same degree-of-freedom. (In the literature, 
this is often called the “drive point” as distinguished from the “transfer” apparent mass.) The 
accelerations at other boundary degrees-of-freedom should be constrained to be zero if one is 
dealing with a multiple drive point problem [8], but herein, only a single drive point is usually 
of interest, and this consideration is ignored. The apparent mass is generally a complex 
quantity, with magnitude and phase, but herein the term apparent mass will often be used in 
referring to only the magnitude. The apparent mass can vary greatly with frequency, as one 
passes through resonances. Therefore the apparent mass reflects the stiffness and damping 
characteristics of a structure, as well as the mass characteristics. 


The closed form solution for the apparent mass of a rod excited at one end and free at the other 
end is given by [33]: 


tan (7toV2co,) - i tanh (7t^to/2to,) 

= M (co) = (ip, c/co) (1+iQ , (6) 

A ( w ) 1 +i tan (tc(o/2co,) tanh (7t^©/2(0,) 

where: p, is the mass per unit length, c is the speed of longitudinal waves (EA/ p, ) l/2 where E is 
Young s modulus and A is cross-section area, £ is the critical damping ratio, and co, is the 
fundamental frequency nc/2L with L the rod length. This result is plotted as the solid line in 
Fig. 4 for a critical damping ratio £ of 2.5%. 


3.3.2 Effective Mass 


Another mass-like quantity of great significance in structural analysis and for impedance 
simulation is the effective mass [15,34], A formal definition of the effective mass, which 
encompasses multiple degrees-of-freedom and off-diagonal terms, as well as equations which 
enable the effective mass to be calculated with NASTRAN, is given in the next section. 

For the beam driven at one end, the apparent mass may also be expressed as a modal expansion 
involving the effective modal masses, m n , and the single-degree-of-freedom frequency 
response factors: 


d+120 

M (CO) = M 0 - I n m n , n = 1,2,3 (7) 

{[(1 - (co /oo n ) 2 ] + i 20 
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Revised 9/27/00 
Equation 7 only 


3.3.1 Apparent Mass 

In this monograph structural impedance will be characterized as “apparent mass”, which is the 
preferred name for the frequency response function (FRF) consisting of the ratio of reaction force 
to prescribed acceleration [32], (Apparent mass symbols will be underscored in this monograph 
to distinguish them from other mass quantities.) The force and prescribed acceleration in the 
apparent mass usually refer to the same degree-of-ffeedom. (In the literature, this is often called 
the “drive point” as distinguished from the “transfer” apparent mass.) The accelerations at other 
boundary degrees-of-freedom should be constrained to be zero if one is dealing with a multiple 
drive point problem [8], but herein, only a single drive point is usually of interest, and this 
consideration is ignored. The apparent mass is generally a complex quantity, with magnitude and 
phase, but herein the term apparent mass will often be used in referring to only the magnitude. 
The apparent mass can vary greatly with frequency, as one passes through resonances. Therefore 
the apparent mass reflects the stiffness and damping characteristics of a structure, as well as the 
mass characteristics. 

The closed form solution for the apparent mass of a rod excited at one end and free at the other 
end is given by [33]: 

F(co) tan (tko/2oJ| ) - i tanh (7t£io/2o>i) 

= M (cu) = (ipi c/to) ( 1 +iQ , (6) 

A(oo) 1 +i tan (mo/2t0|) tanh (7 X^co/2oj j ) 

where: pi is the mass per unit length, c is the speed of longitudinal waves (EA/ pi )' ~ where E is 
Young’s modulus and A is cross-section area, £ is the critical damping ratio, and tOi is the 
fundamental frequency 7tc/2L with L the rod length. This result is plotted as the solid line in Fig. 

4 for a critical damping ratio g of 2.5%. 

3.3.2 Effective Mass 

Another mass-like quantity of great significance in structural analysis and for impedance 
simulation is the “effective mass” [15,34], A formal definition of the effective mass, which 
encompasses multiple degrees-of-freedom and off-diagonal terms, as well as equations which 
enable the effective mass to be calculated with NASTRAN, is given in the next section. 

For the beam driven at one end, the apparent mass may also be expressed as a modal expansion 
involving the effective modal masses, m n , and the single-degree-of-freedom frequency response 
factors: 


(l+i20 

M.(w) = I n m n ,n= 1,2,3 (7) 

{[(1 -((u/co n ) 2 ] + i20 
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Apparent Mass Magnitude, Q=20 
Asymptotic Mass (Q=1) 

Modal Mass 
Residual Mass 


Ratio of Frequency to Fundamental (f/fo) 

FIGURE 4. Apparent Mass, Asymptotic Mass, Modal Mass, and Residual Mass of 

Longitudinally Vibrating Rod, Excited at One End and Free at Other End 
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where the structural form of damping has been assumed, M o is the total mass 
the effective mass of the nth mode which is given by [34]: 


p, L, and m n 


is 


m „ = 8M 0 /[7r(2n-l) 2 ], 


( 8 ) 


where (0 is the natural frequency of the nth mode which is equal to to. (2n -1). Equation 7 mav 
be viewed as a definition of the drive point effective mass. The sum of the effective masses 
over all the modes is the total mass, which may be verified by summing the rod effective 
masses, given by Eq. 8, over all n [34], 


3.3.3 Residual Mass 


A nother mass quantity closely related to effective mass is the “residual mass”, which is defined 

th^ 6 t0 f td t ma f SS mmUS th ? effec ' tlve mass of th e modes which have natural frequencies below 
the excitation frequency. Thus the residual mass of the Nth mode is: 


M (N) = M - 1 , M m 

v 7 o n = 1 to N 11 1 n ‘ 


(9) 


The residual mass may be interpreted as the fraction of the total mass which moves with the 
input acceleration, like a rigid body. A more precise definition of the residual mass concept, 
which encompasses multiple degrees-of-freedom and off-diagonal terms, is given in the next 

It follows from the definition of residual mass, and the fact that the sum of the effective masses 
n; the total mass, that the residual mass decreases monotonically to zero as frequency increases. 
This is the mechanical analogy of Foster’s Theorem for electrical circuits [3 1 ] Herein the 
residual mass is generally indicated with a upper case M. The effective modal mass is the 
negative change in the residual mass, at the resonance frequencies. 

are\1s^ show^n Fig 1 4 fesidual masses of the first five modes of the rod excited longitudinally 

The dashed curve in Fig. 4 is the critically damped apparent mass which may be obtained from 

Eq. 6 by setting the critical damping ratio £ equal to unity. The critically damped apparent mass 

s S vSm” ! ^r e “nl h v m m e t e ^ n funCtlon in electricaI circuit analyses [31] or alternately the “infinite 
system or asymptotic apparent mass. 

The name “infinite system” derives from the second method of calculating this function which 
mi y NJ OI I Slde fv! n f g t sei ™' inf,mte system, e.g. the first factor on the right-hand-side of Eq. 6 
fin ' n TK e , infm,te s y st e m value of apparent mass of the rod decreases as one over 
trequency. This is also true for the apparent mass of a plate vibrating in bending. 

The name asympotic derives from the third method of calculating this function which is to take 
a geometric average of the apparent mass FRF over frequency, so that there is equal area above 

rrilf 1 0W H C T e 0n a log ' log . plot [35 ’ 36 J- (° bserve the equal area characteristic of the 

imnortnninih^H api ] arent mass in Fl S‘ 4 ) The asymptotic form of the apparent mass is very 
important to the development herein, because it will be used to represent experimental apparent 

mass data measured either in tap tests or shaker tests. Notice in Fig. 4 that the asymptotic 
apparent mass is a generous envelope of the residual mass. The asymptotic apparent mass 
which must include stiffness as well as mass contributions, is approximately equal to 2 I/2 times 
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,he residual mass a, the natural frequencies in Fig. 4. Herein the asymptotic apparent mass will 
be used as an approximation to the residual mass. 


3.3.4 FEM Calculation of Effective Mass 

Subdividing the displacement vector into unrestrarned absolute drsplacemen.s u, and presented 

absolute displacements, u p the equilibrium Equation is [15, 34] . 

[ m FF I m FP ] d 2 u F /dt 2 ^ [ k FF I k FP ] u F ^ ^ |f ^ ^ (10) 

[ m PF I m PP ] d 2 u P /dr [k PF lk PP ] u p f P 

L 4>n • 4>p 1 U n For correct printing of Eqs.10-12 

Let: {u} = (j) U = { — }’ see Eqs. B1-B3 in Appendix C 

[ 0 I I PP ] U P 

generalized prescribed absolute displacement. Substituting and pre-multiplying by <D yields. 

[M nn IM np 1 d 2 U N /dt 2 ^ j« 2 _N M NNj _ 0 ]_ U N F. } = { _ _ K (12 ) 
[M t np IM pp ] d : U P /dt 2 [0 I 0] Up Fp 

t w (13) 

where: M NN - (j) N m FF sVi 

M N p = (t)N rm FF ( l ) P + < J ) N Tm FP I PP T (15) 

M PP = I PP m pp I P p + lpp m PF <t> P + W m FP !pp + ^ < * )p ( 1 6) 

F p — lpp fp • 

For: d 2 U P /dt 2 = U P = F f = 0, d 2 U n /dt 2 = - co 2 U n , and for U n = 1 : 

M nP T = - F P / a>„\ ° 7) 

where n indicates a single mode. ta"" - ““ 

elastic-rigid coupling or the modal participa SPCFORCE at that point in a 

restrained at a single point, the reaction (F p ) in Eq. 17 is the 

NASTRAN modal analysis, 
the effective mass: 

, (18) 
M„ T M„ ' M, , 

from the current M PP ", which is the res, dual mass after excluding .he mass associated with the 
already processed n modes. 
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P.'o the prescribed acceleration 

thC rea ? tlon is desired in the same direction as the excitation and the effective mass 
for the common direction is a diagonal of the M T M 1 M m o, Tlf Ct f ™ ass 

m P p 0 the residual mass after processing a complete set of modes fs a 6x6 null matrix Ifm 
and m pp are no, equal ,o zero. ,he vaine of M p /af,er pr«essin g e a c^ptae “e, of modes 

m pp m Fp 0N M NN (J) n m FP , which must be positive definite. 
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4.0 Force Limits 

There a. vrr.ua, ly no high, data , 

structure and test item interfaces. Curr y desree-of freedom or other analytical models 

using one of three methods: 1) mass data for the 

of the coupled source/load to g e t d ■ semi -empirical method based on system test 

mounting structure and test item, 2) estimated u g ic J esi criteria w hich may be based 

data and heuristic arguments, or 3) taken f 9 . In t f, e f irst two methods, which are 

on coupled loads analysis or a simple mas rP pime the force specification is based on and 

usual,/ applicable in the random a “!on : s q “S 8 o rXn^co3sm (or error) in .he 
proportional to the conventional accelerat P { * j the third me thod, which is 

& S “ .e,s, .he force ,inu. ,s 
derived independently from the acceleration specification. 


4.1 Coupled System Methods 

The basic approach to calculating force limits with a coupled system model involves four steps: 

1 Develooment of a parametric model of the coupled source and load system 
‘ STght be anFEM or a multiple degree-of-freedom modal model, 

2. Identification of the model parameters using measured apparent mass or FEM 
modal frequency and effective mass information, 

1 Solution of the coupled system problem to obtain the ratio of the force frequency 
envelop ?he Sera, ion envelope a, the source/load mterface, and 

4, Multiplication of the ratio of envelopes by the acceleration specification to obtain 
an analogous force specification. 

“Sens implemented for'two^E^fic cates' S"- Simple and a 
complex two-degree-of-freedom system (TDFS). 

spectral density S AA as [27]: 


(to) = IM 2 (o)I 2 S aa (to). 


(19) 


feme and acc'eleratton applies a, each frequency. 
It can be shown [37] that, for white noise 

acceleration spectral peak at that natural frequency. 
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asisgs s=s?s=- 

™ S . model sh °-" i" Rg. 2 when >he 
e q ua, ,o the product of theC ^n,“ *, 

4.1.1 Simple TDFS Method 

^"ourc™' nd So C a*t!na,ors p 7 f ^ '" Rg ' 2 With different ™ sses 

therefore the maximum interface feme rienr when I F ' ll,e 1 m “ ,mum response of the load and 

sssksss? [38) - 

(Co/coj- = l+(m 2 /m,)/2 ± [(m^/m, )+(m 2 /m, ) 2 /4)] 05 . ^0) 

rs r tral C * U,attd as 19 from the 


s ff/(S aa rn 2 2 ) = [ 1 +( 0 )/ 0 ) o ) 2 /Q 2 2 ] /{fl-(co/(o 0 ) 2 ] 2 +(qV 0 ) o ) 2 /Q 2 2 }. 


( 21 ) 


S32at&'S=S'*S 

frequency from Eq 20 is plotted aaaincf th,* r 1 y evaluatin g 21 at the lower resonance 
Fig 5. q ’ ’ P agamst the ratl ° of load to SOUf ce mass for three values of Q, in 

SsSSSSSm^sa- 

from the proplrtk^ 2 be deduced 

« is to identify the oscillToTmL^m 
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Ratio of Load Mass to Source Mass (M^,) 

FIGURE 5. Random Vibration Force Specification Calculated From Simple TDFS 


Defining the oscillator masses as residual masses instead of as modal masses, over estimates both 
the load and source masses which over estimates the interface forces calculated using Fig. 5. I his 
approach is conservative for testing, but it is not very conceptually satisfying. The Fig. 2 model 
with only one mass for the source and another for the load, is basically deficient in that it can not 
represent the force contributions of both resonant and non-resonant structural modes. This leads to 
the consideration of a more complex TDFS model, which for some configurations predicts even 
higher force limits than the simple TDFS model. 


4.1.2 Complex TDFS Method 

Here the force limit is calculated for a more complex TDFS model in which the source and load 
each have two masses to represent both the residual and modal masses of a continuous system 
1271. As in the case of the simple TDFS, it is assumed that the acceleration specification correctly 
envelopes the higher of the two acceleration peaks of the coupled source and load system. 
However, for the more complex TDFS, the calculation of the force limit requires evaluating the 
relative sizes of the acceleration peaks at the two coupled system resonance frequencies (i.e. 
determining the mode shapes of the coupled system), which requires some specific assumptions 
about how the system is excited. Calculation of the force limit for this system also necessitates a 
tuning analysis, in which the maximum force is calculated for different ratios of the load an 
source uncoupled resonance frequencies. The complexity of the model requires that the results be 
presented in tabular form for different ratios of modal to residual mass, for both the source and the 

load. 
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Fig. 6a shows a model of a source or a load in which each mode may be represented as a collection 
of single-degree-of-freedom systems attached in parallel to a common interface. (This type of 
model is sometimes called an asparagus patch” model [34].) The masses of the single-degree-of 
freedom systems in this parallel oscillator model are defined by the effective masses of the modes 
of the distributed system. 

When this parallel oscillator model is excited at the interface with a frequency near the resonance 

frequency 0 ) n of the nth mode, the model may be simplified to that in Fig. 6b, where m n is the 
modal mass of the nth mode and M n is the residual mass, i.e. the sum of the masses of all modes 
with resonance frequencies above the excitation frequency. Finally, Fig. 6c shows a coupled 
system with a residual and modal mass model of both the source and the load. The ratio of modal 

to residual mass is 0C]=m|/M| for the source and 0 t 2 = ni 2 /M 2 for the load; the ratio of load and 

source uncoupled resonance frequencies is and the ratio of load and source residual 

masses is (i=M 2 /M,. 

The undamped resonance frequencies of the coupled system in Fig. 6c are solutions of: 


(l-B, 2 )(l-6 2 2 )+ ot,( 1 -B 2 2 ) + |J.( 1 -6, 2 )( 1 -B 2 2 )+^ux 2 ( 1 -B , 2 ) = 0, 

with 8,={fl/C0|, 6 2 =co/co 2 , to,=(k ] /m 1 ) 05 , and (o 2 =(k 2 /m 2 ) 05 . 


( 22 ) 


Using Q to eliminate 6,, the two undamped resonance frequencies of the coupled system are found 
from the quadratic equation solution: 


B 2 2 = -B/2 ± (B 2 -4C) 05 /2, 


(23) 


where: 


B - -[(l+fi+a,)/Q 2 +(l+^+fia 2 )]/(l+|j.) and C = (l+M.+a,+fia 2 )/[(l+|Li)Q 2 ]. 

The interface force spectral density, normalized by the acceleration spectral density and the load 
residual mass squared, calculated as in Eq. 19 from the magnitude squared of the load dynamic 
mass, is: 


IM//M 2 2 = { [( 1 -B 2 2 )+a 2 ] 2 +8 2 2 ( 1 +a 2 ) 2 /Q 2 2 }/ [(l-B 2 2 ) 2 -t-6 2 2 /Q 2 2 ] . 


(24) 


Combining Eqs. 23 and 24, yields the normalized force spectral density at each of the two coupled 
system resonance frequencies. 

The desired result is the ratio of the larger of the two force spectral density peaks to the larger of 
the two acceleration spectral density peaks, the former being the desired force limit and the latter 
corresponding to the acceleration specification. Unfortunately, the peak acceleration and peak force 
do not necessarily occur at the same frequency, e.g. the peak acceleration may be at the higher of 
the two coupled system resonance frequencies while the peak force is at the lower of the two 
frequencies. This, in fact, occurs when the uncoupled resonance frequencies of the load and source 
are approximately equal, that is for D. near unity. 
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ia/a Io i 2 = im,/(M i +M 2 )I 


(25) 


«££ &£££ %£SZ£SZfc£Sf££* M ; wh,ch r 

source modal mass m,, the magnitude squared of the free acceleration A lo is: e ing ° n the 


IA lo /(F e /m,)l 2 = B/O+BjVq, 2 ) / { [ ( l-8 1 2 )(l-8 1 2 /a 1 )-l] 2 +B l 6 (l + l/a,) 2 /Q l 


; }. (26) 


c^r m - yidds the & 

corresponding to the" mailed acceleration JSk'bythf ratoof tte ^ifenothVr 0 "™” f |' eqU ? nc >' 

S3® d S“Slt„ o r X ncy ? y f y 

the greater acceleration spectral density peak. f H & * f ° KC s P ectral densit y P ea k to 

The final step in the derivation of the force limit is to vary the ratio O = rn / m n fth 

maximum tuned values of the force spectral density Vonri^l inc '[ e h me ^ ts in the frequency ratio. The 

and the maximum value of the acceleration spectraf densatWwhth^ ^ loadresiduaJ mass squared 
specification) are listed in Tables 1 2 anH th<* >ty (which is equivalent to the acceleration 

50, 20, and 5 fact< f bo,h equal 10 

deemed necessary ) J ' K h Q may be com P ut ed from Eqs. 22 -26, if 

The normalized force speclral density in Tables 1-3 is unity when a, = 0 (see Eq. 24) and should 
e interpolated for 0. 1 25 < a, < 0. It is suggested that the force spectral density for a < 0 25 be 

"“Cuu£7 Pt °' ,C ValUe f ° r = °- WhiCh » a ve, T small mass moving a, 
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Table 1. Force Limit Spectrum for Complex TDFS with Q=50 




Ratio of modal 
to residual mass 
mt/MI , m2/M2 


Residual mass ratio 
O.CM 0.03 0.1 


8.0, 

8.0 

8.0, 

4.0 

8.0, 

2.0 

8.0, 

1.0 

8.0, 

0.5 

8.0, 

0.25 

8.0, C 

). 1 25 

8.0, 

0.0 

4.0, 

8.0 

4.0, 

4.0 

4.0, 

2.0 

4.0, 

1.0 

4.0, 

0.5 

4.0, 

0.25 

4.0, 0.125 

4.0, 

0.0 

2.0, 

8.0 

2.0, 

4.0 

2.0, 

2.0 

2.0, 

1.0 

2.0, 

0.5 

2.0, 

0.25 

2.0, 

0.125 

2.0, 

0.0 

1.0, 

8.0 

1.0, 

4.0 

1.0, 

2.0 

1.0, 

1.0 

1.0, 

0.5 

1.0, 

0.25 

1.0, 

0.125 

1.0, 

0.0 

0.5, 

8.0 

0.5, 

4.0 

0.5, 

2.0 

0.5, 

1.0 

0.5, 

0.5 

0.5, 

0.25 

0.5, 

0.125 


0.5, 0.0 

0.25, 8.0 
0.25, 4.0 
0.25, 2.0 
0.25, 1.0 
0.25, 0.5 
0.25, 0.25 
0.25, 0.125 
0.25, 0.0 

0.125, 8.0 
0.125, 4.0 
0.125, 2.0 
0.125, 1.0 
0.125, 0.5 
0.125, 0.25 
0.125, 0.125 
0.125, 0.0 

0 . 0 , 8.0 
0.0, 4.0 
0.0, 2.0 
0 . 0 , 1.0 
0.0, 0.5 
0.0, 0.25 
0.0, 0.125 
0.0, 0.0 


949 

950 
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237 

237 
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59 

59 

59 

1 5 

1 5 

1 5 


4 4 4 
1 1 1 
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870 

221 
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220 

55 

55 

55 

1 4 

1 4 

1 4 

3 

3 

3 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1640 

1521 

1286 

420 

404 

364 

1 06 

105 

99 

27 

27 

26 

7 

7 

7 

2 

2 

2 

1 

1 

1 

1 

1 

1 

16554 

6508 

2921 

7333 

2965 

1200 

3080 

1 345 

502 

1 1 89 

592 

229 

41 5 

245 

106 

1 32 

94 

48 

39 

33 

2 1 

1 

1 

1 

24199 

9798 

3726 

10238 

4417 

1672 

4046 

1927 

747 

1 454 

804 

335 

472 

31 1 

1 48 

141 

1 1 0 

63 

40 

36 

26 

1 

1 

1 

3391 0 

13269 

4455 

14189 

6185 

2155 

5342 

2736 

1043 

1 764 

1111 

492 

529 

396 

219 

149 

128 

85 

4 1 

39 

3 1 

1 

1 

1 

48146 

18637 

6361 

19122 

8823 

2885 

6642 

3788 

1 454 

2045 

1434 

684 

574 

477 

291 

155 

1 42 

1 1 0 

4 1 

4 1 

36 

1 

1 

1 


134767 

66196 

13561 

37885 

28769 

8669 
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8998 

5015 

2484 
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625 
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40 

40 

40 

1 

1 
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59 

60 
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4 

4 
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1 

1 
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1 

1 

1 

1 

860 

916 
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57 
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4 

4 
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1 
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887 
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89 

1 42 

62 

40 

66 

30 

1 8 

3 1 

1 5 

9 

1 5 

8 

5 
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393 
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405 
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96 
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80 

45 

104 

45 

23 

47 

22 
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20 

1 1 

6 

1 

1 

1 

241 1 
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41 1 
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508 
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1 05 
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1 39 
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66 
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27 
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1 0 

1 

1 

1 

2836 

1 1 36 
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1 827 
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258 
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1 1 0 
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1 87 

54 
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1 1 1 

35 
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69 

25 

38 

30 

1 4 
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22 
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8 

8 
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5 

5 

5 

1 

1 

1 

936 
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1 194 
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244 
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74 

67 

66 

30 

24 

22 

1 4 

1 3 

1 1 

9 

7 

7 

6 

5 

5 

1 

1 

1 

917 

1110 

119 1 

233 

245 

245 

69 

68 

66 

30 

25 

22 

1 7 

1 2 

1 1 

1 0 

8 

7 

9 

5 

4 

1 

1 
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Table 2. Force Limit Spectrum for Complex TDFS with Q=20 

(Normalized By Load Residual Mass Squared and Acceleration Spectru m) 

Ratio of modal Residual mass ratio , M2/M1 

to residual mass 0.001 6.003 6.01 oHo3 o"l O i 

ml /Ml . m2/M2 


8.0, 8.0 

932 

933 

936 

948 

8.0, 4.0 

233 

233 

233 

235 

8.0, 2.0 

58 

58 

58 

5 8 

8.0, 1.0 

1 5 

1 5 

1 5 

1 5 

8.0, 0.5 

4 

4 

4 

4 

8.0, 0.25 

1 

1 

1 

1 

8.0, 0.125 

1 

1 

1 

1 

00 

o 

o 

o 

1 

1 

1 

1 

4.0, 8.0 

871 

867 

858 

849 

4.0, 4.0 

21 8 

218 

217 

21 6 

4.0, 2.0 

55 

55 

55 

55 

4.0, 1.0 

1 4 

1 4 

1 4 

1 4 

4.0, 0.5 

3 

3 

4 

4 

4.0, 0.25 

1 

1 

1 

1 

4.0, 0.125 

1 

1 

1 

1 

o 

b 

o 

1 

1 

1 

1 

2.0, 8.0 

1 586 

1 478 

1260 

1061 

2.0, 4.0 

406 

391 

355 

305 

2.0, 2.0 

1 03 

101 

97 

88 

2.0, 1.0 

26 

26 

26 

25 

2.0, 0.5 

7 

7 

7 

7 

2.0, 0.25 

2 

2 

2 

2 

2.0, 0.125 

1 

1 

1 

1 

2.0, 0.0 

1 

1 

1 

1 

1.0, 8.0 

1 1041 

5731 

2714 

1486 

1.0, 4.0 

3869 

2206 

1 1 05 

567 

1.0, 2.0 

1228 

826 

432 

226 

1.0, 1.0 

359 

283 

166 

100 

1.0, 0.5 

1 00 

89 

63 

42 

1.0, 0.25 

28 

27 

23 

1 7 

1.0, 0.125 

8 

8 

8 

7 

1.0, 0.0 

1 

1 

1 

1 

0.5, 8.0 

1 3889 

7720 

3501 

1 726 

0.5, 4.0 

4516 

2895 

1417 

695 

0.5, 2.0 

1 346 

1 003 

561 

283 

0.5, 1.0 

377 

31 9 

21 1 

1 1 7 

0.5, 0.5 

1 02 

95 

74 

48 

0.5, 0.25 

28 

27 

25 

1 9 

0.5, 0.125 

8 

8 

8 

8 

0.5, 0.0 

1 

1 

1 

1 

0.25, 8.0 

1 7378 

9978 

4092 

1 944 

0.25, 4.0 

5194 

3725 

1 805 

812 

0.25, 2.0 

1 455 

1 205 

741 

359 

0.25, 1.0 

391 

354 

269 

160 

0.25, 0.5 

1 03 

99 

86 

63 

0.25, 0.25 

28 

28 

27 

23 

0.25, 0.125 

8 

8 

8 

8 

0.25, 0.0 

1 

1 

1 

1 

0.125, 8.0 

19966 

12425 

5389 

2331 

0.125, 4.0 

5748 

4417 

2241 

1080 

0.125, 2.0 

1 533 

1 368 

901 

429 

0.125, 1.0 

400 

380 

312 

1 92 

0.125, 0.5 

1 04 

102 

95 

75 

0.125, 0.25 

27 

28 

27 

26 

0.125, 0.125 

8 

8 

8 

8 

0.125, 0.0 

1 

1 

1 

1 

0.0, 8.0 

251 14 

21284 

10111 

2700 

0.0, 4.0 

6394 

6108 

4156 

1 560 

0.0, 2.0 

1608 

1 590 

1409 

757 

0.0, 1.0 

404 

403 

390 

310 

0.0, 0.5 

1 02 

1 02 

101 

95 

0.0, 0.25 

27 

27 

27 

26 

0.0, 0.125 

8 

8 

8 

7 

o 

o 

o 

b 

1 

1 

1 

1 
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1 180 

1 240 

1234 

1 238 

239 

256 

294 

265 

250 

59 

60 

68 

73 

68 

1 5 

1 5 

1 7 

23 

22 

4 

4 

4 

7 

6 

1 

1 

1 

2 

5 

1 

1 

1 

1 

3 

1 

1 

1 

1 

1 

904 

1042 

1 067 

1110 

1229 

220 

250 

254 

250 

252 

56 

61 

72 

68 

67 

1 4 

1 6 

2 1 

23 

22 

4 

4 

6 

1 0 

1 0 

1 

1 

2 

5 

5 

1 

1 

1 

3 

3 

1 

1 

1 

1 

1 

990 

946 

982 

1099 

1 201 

272 

259 

238 

236 

254 

79 

82 

70 

65 

62 

24 

25 

25 

23 

22 

7 

9 

1 0 
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1 0 

2 

3 

5 

5 

6 

1 

1 

3 

3 

4 

1 

1 

1 

1 

1 

967 
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984 

1095 

1 181 

332 

247 

233 

238 

248 

125 

83 

71 

66 

64 

50 

34 

26 

23 

23 

24 

1 5 

1 2 

1 1 

1 1 

1 1 

8 

6 

6 

6 

5 

5 

4 

4 

4 

1 

1 

1 

1 

1 

1023 

880 

974 

1093 

1 1 71 
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247 

225 
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244 

1 36 

89 

70 

64 

65 

59 

39 

27 

24 

22 

27 

1 7 

1 2 

1 1 

1 0 

1 3 

8 

7 

6 

6 

6 

5 

4 

4 

4 

1 

1 

1 

1 

1 

1017 

833 

936 

1092 

1 166 

380 

249 

225 

241 

242 

1 73 

93 

71 

66 

65 

74 

43 

28 

23 

22 

38 

22 

1 4 

1 2 

1 1 

1 6 

1 0 

8 

7 

7 

7 

5 

5 

4 

4 

1 

1 

1 

1 

1 

1048 

849 

918 

1092 

1 1 63 

400 

266 

228 

242 

241 

1 84 

1 02 

72 

66 

65 

91 

45 

29 

24 

22 

42 

24 

1 4 

1 2 

1 1 

20 

1 3 

8 

7 

6 
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1 

1 

1 

1 

1 

1125 
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900 
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240 
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1 09 

68 

67 

66 

1 48 

52 

30 

25 

22 

60 

30 
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1 2 

1 1 

22 

1 7 

9 

7 

6 

7 

6 

6 

5 

4 

1 

1 

1 

1 
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Ratio of modal 
to residual mass 
ml/MI . m2/M2 

8 . 0 , 8.0 
8.0, 4.0 
8 . 0 , 2.0 
8 . 0 , 1.0 
8.0, 0.5 

8.0, 0.25 

8.0, 0.125 
8.0, 0.0 

4.0, 8.0 
4.0, 4.0 
4.0, 2.0 
4.0, 1.0 
4.0, 0.5 
4.0, 0.25 

4.0, 0.125 

4.0, 0.0 

2 . 0 , 8.0 
2.0, 4.0 
2.0, 2.0 
2.0, 1.0 

2.0, 0.5 

2.0, 0.25 

2.0, 0.125 

2 . 0 , 0.0 

1.0, 8.0 
1.0, 4.0 
1 . 0 , 2.0 
1.0, 1.0 
1.0, 0.5 

1.0, 0.25 

1.0, 0.125 

1.0, 0.0 

0.5, 8.0 
0.5, 4.0 
0.5, 2.0 
0.5, 1.0 
0.5, 0.5 
0.5, 0.25 
0.5, 0.125 
0.5, 0.0 


Table 3. Force Limit Spectrum for Complex TDFS with Q=5 

(Norm alized By Load Residual Mass Squared and Acceleration Spectrum) 

' Residual mass ratio , M2/M1 

0.001 0.003 0.01 0.03 (TT 0*3 i ~ 
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1 77 

1 78 

1 78 
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46 

46 

46 

46 

1 2 

1 2 

1 2 

1 2 

4 

4 

4 

4 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

687 
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685 

689 

1 74 

1 74 

1 74 

1 75 

45 

45 

45 

45 

1 2 

1 2 

1 2 

1 2 

3 

3 

3 

4 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 006 

983 

927 

860 

256 

254 

247 

235 

66 

66 

65 

64 

1 8 

1 8 

1 8 

1 8 

5 

5 

5 

5 

2 

2 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 700 

1618 

1389 

1048 

434 

430 

397 

324 

1 1 3 

1 1 3 

1 1 1 

102 

30 

31 

31 

31 

9 

9 

1 0 

1 0 

3 

3 

4 

4 

2 

2 

2 

2 

1 

1 

1 

1 

1 703 

1 655 

1456 

1 1 47 

433 

432 

409 

350 
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1 1 3 

1 12 
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31 

3 1 

3 1 
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1 0 
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1 0 

4 

4 

4 

4 

2 

2 

2 

2 

1 

1 

1 
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752 

866 
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876 

183 

1 97 

221 

21 2 

46 

48 

56 
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1 2 
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1 9 

4 

4 

4 
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1 

1 

1 
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1 

1 
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1 

1 

1 

1 

1 
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827 
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1 82 
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203 
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52 

59 
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4 

4 
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1 

1 

2 

4 

1 

1 

1 

2 

1 

1 

1 

1 
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64 

64 

57 

56 
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2 1 

20 
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6 

7 

8 

8 

2 

3 

4 

4 

1 

1 

2 

2 

1 

1 

1 

1 

787 
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835 
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80 

66 

58 

55 

28 

24 

20 

20 

1 0 

9 

8 

8 

4 

4 
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4 
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2 

2 
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2 

2 
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0.25, 

8.0 

1698 

0.25, 

4.0 

434 

0.25, 

2.0 

1 1 2 

0.25, 

1.0 

31 

0.25, 

0.5 

1 0 

0.25, 

0.25 

4 

0.25, 0.125 

2 

0.25, 

0.0 
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To use the force limits in Tables 1-3, both the residual and modal masses of the source and load 
must be known as a function of frequency, either from a test, from an FEM, or from both. FEM 
analyses provide both the modal and residual effective masses (see Section 3.3.4). If shaker or tap 
tests are used to measure the effective masses, the smoothed FRF of the magnitude of the ratio of 
force to acceleration should be taken as the effective residual mass, as discussed in Section 3.3.3. 


4.2 Semi-empirical Method of Predicting Force Limits 


4.2.1 Rationale 


The semi-empirical approach to deriving force limits is based on the extrapolation of interface force 
data for similar mounting structure and test items. The following form of semi-empirical force limit 
for sine or transient tests was proposed in 1964 [3]: 

F, = C M 0 A s (27) 

where F, is the amplitude of the force limit, C is a frequency dependent constant which depends on 
the configuration, M 0 is the total mass of the load (test item), and A s is the amplitude of the 
acceleration specification. The form of Eq. 27 appropriate for random vibration tests is: 

S FK = C 2 M 0 2 S AA (28) 

where S n . is the force spectral density and S AA the acceleration spectral density. 

In [3], it is claimed that C seldom exceeds 1.4 in coupled systems of practical interest because of 
the vibration absorber effect. From the preceding analysis of TDFS’s, it is apparent that this claim 
implies something about the ratio of the load and source effective masses of the coupled system. 
For the simple TDFS shown in Fig. 2, the ordinate of Fig. 5 may be interpreted as C 2 in Eq. 28. 
Thus a value of C of 1.4 corresponds to an ordinate value of 1.96 and an abscissa value the ratio 
of load to source mass (M 2 /M,), of 1.5. Although the load and source effective masses’ are often 
comparable in aerospace structures, such is not always the case, and the use of Eqs. 27 or 28 with 
C — 1 .4 might result in undertesting lor lightweight loads mounted on heavy structure. 

A refinement of Eq. 28 follows from inspection of Eq. 19, which is Newton’s second law for 
random vibration. If one takes the envelope of both sides of Eq. 19, the left-hand-side is the 
envelope of the interface force spectrum, which is the sought-after force limit. On the right-hand- 
side results the envelope of the product of the load apparent mass and the interface acceleration 
spectrum. This product may be approximated as the frequency average (the aforementioned 
asymptote) of the magnitude squared of the load apparent mass, times the envelope of the interface 
acceleration spectrum. This refinement of the semi -empirical approach was the subject of several 
important papers in the 1970’s [10,13 & 14], Herein, it will be assumed that the asymptote of the 
magnitude of the load apparent mass is equal to the total mass below the first resonance, and then 
falls off as one over frequency at frequencies above the first resonance, as is the case for the rod 
example in Fig. 4 and for a plate excited in bending. (However, the asymptotic mass of a beam 
excited in bending, falls off as one over the square root of frequency.) Assuming a one-over- 
frequency fall-off of the asymptotic load mass, leads to the following modification of Eq. 28 
applicable to random vibration testing: ^ 
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S FF = C 2 M (1 2 S AA , f < f D 

S FF = C 2 M 0 2 S AA / (f / f„) 2 f>f 0 


( 29 ) 


Some judgment and reference to test data for similar configurations must be considered to choose 
the value of C and the exponent of (f/f 0 ) Eq. 29. Herein are considered force data measured at the 
interface between three equipment items and the Cassini spacecraft Development Test Model 
(DTM) during acoustic tests. 


4.2.2 Validation 

Figure 7 shows a schematic of the Cassini spacecraft which will be launched in October 1997 to 
arrive at Saturn in 2004. Figure 8 shows the DTM spacecraft configured for one of several acoustic 
tests in the JPL reverberant acoustic chamber [39]. The three equipment items which were 
instrumented with tri-axial force transducers between their mounting structure and the spacecraft 
ring-stringer structures are: a dynamic model of a Radioisotope Thermoelectric Generator (RTG) 
three of which provide the spacecraft electrical power, the engineering model Radio Plasma Wave 
Subsystem (RPWS), and the flight Propulsion Module Subsystem Electronic Assembly (PMSEA). 
The RTG weighted approximately 120 lb. and was cantilevered outward from a mounting bracket 
attached at four points to the spacecraft lower equipment module. The RPWS weighed 
approximately 65 lb. and was attached with three trusses to the spacecraft upper equipment 
module. The PMSEA is a large electronic box which also weighed approximately 65 lb. and was 
mounted at four comers to the propulsion module. 

Figure 9 shows the spectra of the total vibration force acting in three directions at the 
PMSEA/spacecraft interface during the DTM protoflight level acoustic test. Figure 10 shows the 
acceleration spectra measured at the PMSEA/spacecraft interface in the acoustic test. Figure 1 1 
shows the flight PMSEA mounted on a shaker for a vertical vibration test. Figure 12 shows the 
magnitude (reduced by 4) of the PMSEA apparent mass measured in a preliminary low-level 
(0.25 G) sine sweep vertical axis vibration test. Notice that the fundamental resonance of the 
PMSEA in the shaker vertical test is at approximately 400 Hz, whereas the fundamental (radial) 
resonance of the PMSEA mounted on the spacecraft DTM (Fig. 9) is approximately 100 Hz. This 
disparity between the fundamental resonance frequencies on the flight structure and on the shaker 
is typical, and it should discourage anyone from thinking that the shaker test is a close replica of 
the field environment. The goal of force control is simply to limit the maximum force on the 
shaker, to that estimated for flight, and it must be recognized that the frequencies at which this 
maximum force is exhibited will be different in the two configurations. 

Also shown in Fig. 9 are the semi-empirical vertical force specification calculated from Eq. 29 and 
the force limit used in the vertical vibration test. In the case of the PMSEA, the vibration test force 
specification was based on an envelope of the DTM acoustic test measurements. (It is very unusual 
to have system acoustic test data available before the instrument vibration test.) The semi-empirical 
force specification is based on the acceleration spectrum in Fig. 10 and a value of C of unity. 
Because the acceleration specification in Fig. 10 drops at 250 Hz, the semi-empirical force limit is 
about 7 dB less than the vibration test force limit at 400 Hz. The semi-empirical force limit is a 
reasonable envelope of the acoustic test force data, particularly in the mid-frequency range where 
structure-borne random vibration is the dominant source. It should be noted that the acoustic test 
data as well as the semi-empirical force limit, which is proportional to the acceleration 
specification, have a 4 dB margin over the predicted flight environment. Fig. 13 shows the 


4-11 



notching, approximately 10 dB, that resulted in the PMSEA vertical random vibration test, when 
the shaker force was limited to the vibration test specification in Fig. 9. 


Figure 14 compares the semi-empirical and vibration test vertical force specifications for the 
Cassini RPWS instrument with the interface force data measured in the DTM acoustic test. Both 
force specifications envelope the interface force data peaks at approximately 45 Hz in the radial 
direction and at 65 Hz in the lateral directions. Figure 15 shows the interface acceleration data 
measured at the RPWS/DTM interface. One of the lateral acceleration measurements greatly 
exceeds the specification. Figure 16 shows the RPWS configured for a lateral vibration test. Figure 
17 shows the magnitude (reduced by 4) of the RPWS apparent mass measured in a preliminary 
low-level (0.25 G) sine sweep vertical axis vibration test. Notice that the fundamental resonance of 
the RPWS in the shaker vertical test is at approximately 250 Hz, whereas the fundamental (radial) 
resonance of the RPWS mounted on the spacecraft DTM (Fig. 14) is approximately 45 Hz. This 
discrepancy explains why the semi-empirical force specification does not roll off until the first 
resonance frequency on the shaker and why the force specifications in Fig. 14 must greatly exceed 
the DTM data in the high frequency regime. Figure 18 shows the notching that resulted in the full- 
level vertical random vibration test of the RPWS, when the test force specification in Fig. 15 was 
utilized. 

Figures 19 and 20 show the Cassini RTG interface forces and accelerations measured at the 
RTG/DTM interface. The acceleration test specification in Fig. 20 is a very accurate envelope of the 
DTM acoustic test data. The first resonance in the radial direction on the spacecraft DTM is at 
approximately 220 Hz. (The radial direction on the spacecraft is along the RTG axis.) The semi- 
empirical force limit in Fig. 19 is flat to 750 Hz because the first RTG resonance in the vertical 
vibration test shown in Fig. 21 was at approximately 750 Hz. The RTG’s were inherited hardware 
which had been previously qualified to an existing test specification. The force specification in Fig. 
19 was used to justify the extension of the previous qualification testing to the Cassini program. 
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FIGURE 7. Cassini Spacecraft Schematic (Showing subject PMSEA, RPWS, and RTG’s) 
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FIGURE 8. Cassini DTM Spacecraft in Acoustic Chamber 
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FIGURE 1 1 . Cassini PMSEA In Vertical Vibration Test 
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FIGURE 12. Force in 1/4 G Sine-sweep Vertical Vibration Test of Cassini PMSEA 



FIGURE 13. Notched Acceleration Input in Vertical Random Vibration Test of PMSEA 
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FIGURE 16. Cassini RPWS In Lateral Vibration Test 
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FIGURE 21. Cassini RTG In Vertical Vibration Test 


4.3 Design Load Verification 

In the low frequency regime of sine and transient vibration testing, the design loads may often be 
utilized for force limiting. In fact, the objective of low frequency vibration testing is usually to 
verify that the structure will survive the loads to which it was designed. With the advent of 
vibration torce measurement and limiting, it is cost effective in some cases to use low frequency 
vibration testing to replace the traditional static testing. In many programs, the structural design 
verification is accomplished by analysis, using higher design margins than would be used for a 
structure which will be verified by test. Use of a low frequency vibration test for design 
verification will reduce the costs of analysis and enable the structure to be designed with less 
margin, which can be translated into cost and weight savings and/or into increased performance. 
As in the case ot high frequency random vibration tests, the test force limit in design verification 
tests usually incorporates some margin over the expected flight load. For example, if the expected 
flight limit is taken as unity, the test maximum load might be 1 .2, and the design load 1.5. 


4.3.1 Quasi-static and Coupled Loads 


Early in the program, the design loads are often given as a “quasi-static” acceleration of the CG of 
the test item. The quasi-static loads typically combine the steady loads associated with the 
acceleration of the rocket at launch with the low frequency transients due to launch vehicle staging 
and low frequency aerodynamic loads associated with gusts, buffeting, etc. The quasi-static loads 
are interpreted by the designer as static loads. Therefore these quasi-static accelerations are 
multiplied by the total mass of the test item to obtain the applied loads used for preliminary design, 
for selecting the fasteners, etc. These quasi-static accelerations are obtained from past experience or 
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from semi-empirical curves developed for different launch vehicles and spacecraft configurations 
[40], 

Later in the program, the low frequency loads used for design ^ 

designs are refined by analysis of coupled finite element models (FEM s) [4 1] The spacecratt 
loads are determinedfrom a coupled spacecraft and launch vehicle model, and later the equipment 
loads are determined from a coupled equipment and launch vehicle model. Obviously the vali i y 
oftc FEM Sels depends on the skill of the modeler and upon the accuracy of the information 
used to develop the model. More sophisticated models are used as the programs mature. For 
example Inpreliminary analyses of a spacecraft, much of the equipment may be modeled as a 
lumped masses at theoquipment CG’s. In this case, the calculated forces at the spacecraft and 
equipment interfaces are only valid at frequencies below the first resonance frequency of 
eauinment Later in the program, the equipment may be represented by complex FEM s. B 
these models are usually valid only in the frequency regime encompassing the first few mode, 
the eauinment in each axis. One objective of coupled loads analyses is to maximize the upper 
frequency limit of the model. The validity of the models is often verified and extended to hig er 
frequencies by refining the models with modal test data. 

4,3.2 Force Transducer Measurement of CG Acceleration 

Although it is relatively easy to predict the CG acceleration, it is difficult or impossible t to measure 
the acceleration of the CG with accelerometers in a vibration test. Sometimes the CG 
inaccessible or there is no physical structure at the CG location on which to mount an 
accelerometer However there is a more serious problem. Only in the case of a rigid o y 
CG a fixed point on the structure. Once deformations and resonances and occur, it is impossible 
measure theCG acceleration with an accelerometer. Unfortunately, attempts to measure the CG 

acceleration with an accelerometer usually overestimate the CG r f P ons \ a ^ ^Toaccele ^ion H 
measurements to the CG criterion will result in an undertest. However, the CG acceleration is 

uniquely determined by dividing an interface force measurement by the total mass of the test item, 
per Newton’s second law. 

The non-fixitv of the CG of a deformable body is demonstrated with an example in Fig- 22, which 
^Imtratestlie^lhrd vibration mode of a three-mass, two-spring vibratory system. The mass value 
oUhe mTddle miss I s twice that of the end masses and the two springs are tdenttcal. The upper 
iteKhta Fig shows the system at rest with the CG elearly located at the center of the huddle 
mass The lower sketch in Fig. 22 shows the system displaced m its third mode, with the middle 
mass moving one unit to the left and the two end masses both moving one unit to the right. (The 
first mode involves rigid body translation, and the second mode involves zero motion of the middle 
mass and the two end^nasses moving an equal amount in opposite directions.) It is a ch^teristi 
rfmodal motion that there are no external forces acting, 

displacement illustrated in Fig. 22 can’t involve motion of the C G- However since ™ e ™ 1C £ 
mass, as well as the end masses, are moving, the CG is not at a fixed point in the system. Clearly, 

one could not attach an accelerometer at the CG position. 

The difficulty of measuring the CG acceleration with an accelerometer is further illustrated with 
data obtained on SSSSwl instrument (Fig. 16) in the Cassini spacecraft DTM acou*£ (Fig- 
Ficmre 23 is a schematic of the Cassini RPWS instrument shown in Fig. 16. In the Cassini 
DTM S spacecraft acoustic test, the RPWS was instrumented with tri-axial force transducers be ween 
the instrument and spacecraft. (See force data in Fig. 14.) In addition to the interface 
accelerometers, them was also a tri-axial accelerometer located approxtmately at the CG °f th 
RPWS instrument, in the spacecraft DTM acoustic test. (See position 23 in Fig. 23.) Fig. 24 
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rvr\jf S the ratlo ° f tbe tota external radial force to the radial CG acceleration for the RPWS in the 
DTM spacecraft acoustic test, and Fig. 25 shows the corresponding ratio for a lateral axis. If the 
n ' a ^‘ a ^ a oc e lorometer actually measured the CG acceleration, the curves in both Figs 24 and 05 

ITth^n 6 flat f rr,^ f freque " cy ^ quaI to the totaI wei § ht of the RPWS, approximately 65 lb. 

1 h n P f CCrat i radia i direct,on F «g 24 shows that the ratio offeree to acceleration falls off 
rapidly above about J60 Hz (Fig. 17 shows that there is an RPWS resonance at approximately 
160 Hz in the RPWS vertical direction which corresponds to the spacecraft radial direction ) In the 
spacecraft lateral direction, Fig . 25 shows that the ratio falls off above 50 Hz, which corresponds 
to a lateral resonance of the RPWS on the spacecraft. That the measured ratios of force to 
accderation are less than the total weight indicates that the measured accelerations are greater than 
the true CG acceleration at the higher frequencies. B 

The examples in Figs. 22-25 demonstrate the problems and dangers of using an accelerometer to 
™ a f!f acc ^ leratlon of the CG in vibration tests. It is for this reason that CG response limiting and 

" tlC deSlgn rS fiC r n ha f been difflCult ’ at best ’ to implement previously in vibration 
tests. However with the advent of piezoelectric, tri -axial force transducers, these measurements 

attrac^ve Strai§ht f0rWard> and ^ USC ° f vlbratlon tests for design verification becomes very 


4.4 Force Limiting Vibration Example-ACE CRIS Instrument 

Figure 26 is a photograph of the ACE spacecraft Cosmic Ray Isotope Spectrometer (CRIS) 
instrument mounted on a shaker for a vertical vibration test. The instrument is mounted on twelve 
forces 13 *° rCe transducers ’ wbicb sta y with the instrument in flight and record the lift-off vibratory 

romntv^nPQ Prea m S a Cet f ? f u alCUlating the force limits three wa y- s; the simple TDFS method, the 
5?7f^h7 DFS i Tnpc ar l d tbe semi empirical method. The spreadsheet is linked to Eqs. 20 and 
cointei THF? e ^^calculation and to Tables 1-3, plus an interpolation routine, for the 
complex TDFS calculation. A value of Q of 50, 20 or 5, corresponding to Tables 1 2 and 3 must 

bLdsTom^O m 2? T?u S J alculatlon - The force limits are calculated in one-third octave 
bands from 40 to 1000 Hz. (This frequency range is typical for an instrument, but might include 

h°.^Jk f ~ quenci 5? “J the case of a spacecraft test.) There is nothing sacred about one-third octave 
bands th ° f ne c t ould cb ° ose one-tenth or alternately, octave bands, since specifications are relatively 
smooth functions of frequency. (It should be noted that the width of the notching is set by the 
shaker controller analysis bandwidth, typically 5 Hz for random vibration tests, not by the 
specification bandwidth.) The input acceleration specification must be entered for each one-third 

mem^ ba fTh Reca ^ at th f l° rce hnaUS P r °P ortional to the acceleration specification for all three 
methods.) The remainder of the inputs in the spreadsheet deal with the structural impedance 
characteristics of the load (the test item) and of the source (the mounting structure) For both the 
load and source, one must enter the residual weight and the number of modes in each one-third 
octave band. The spreadsheet calculates modal weight, and then the force limits are automatically 
calculated using the three aforementioned methods. The residual weight information may be Y 
determined from test, FEM, or a combination of these. In this example the information was 
determined from test data. 

Figure 28 shows the magnitude of the apparent mass of an ACE spacecraft honeycomb panel 
measured in a tap test at one of the CRIS instrument mounting locations. (This type of test involves 
tapping at the selected point on the panel near an accelerometer with a small hammer which 
incorporates a force transducer, and using a two channel frequency analyzer to compute the 
magnitude of the apparent mass.) The apparent mass data are smoothed in frequency to compute 
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the asymptotic value, which is taken as the residual mass, as discussed in section 3.3.3. The 
asymptotic value in Fig. 28 rolls off as one over frequency squared above 50 Hz, which is 
characteristic of a spring. Below 50 Hz, the apparent mass looks like a mass, but it is known that 
the coherence fell off below 50 Hz, so the data below 50 Hz are suspect. To obtain the residual 
mass of the source for the spreadsheet, the data in Fig. 28 was multiplied by four to approximately 
account for the multiple (twelve) mounting points. (At low frequencies, the residual mass should 
approach the total mass or stiffness.) In the Fig. 27 spreadsheet, it is assumed that there is one 
source mode in every one-third octave band, so the decrease in the residual mass in each band 
becomes the average modal mass in that band. (When FEM information is used, the number of 
significant modes in each one-third octave band are counted.) 

Figure 29 shows the magnitude of the apparent mass of the CRIS instrument measured on the 
shaker in Fig. 26. The apparent mass has the characteristics of a vibrating plate or rod driven at a 
point, i.e. it is equal to the total mass below the first resonance frequency and then rolls off like 
one over frequency. (See section 3.3.3.) The asymptotic values in Fig. 29 are multiplied by four to 
scale to one G of input, and entered as the load residual mass in the spreadsheet. (Notice that the 
roll-off of the asymptotic mass is started one-third octave band below the resonance frequency in 
Fig. 29 to generate some modal mass at the first resonance frequency.) It is also assumed that there 
is one mode of the load in every one-third octave band above the first resonance frequency. 

The plot in Fig. 27 shows the force limit calculated with each of the three methods, using a 
constant C = 1.5 in Eq. 29 for the semi-empirical method. In this case, the semi-empirical method 
approximately splits the difference between the two TDFS methods, with the simple TDFS method 
giving higher limits below the first resonance and the complex TDFS method giving higher limits 
above the first resonance. 

Figure 30 shows the total vertical force limited to a force specification in the protoflight vertical 
random vibration test of the CRIS instrument shown in Fig. 26. The force specification actually 
used in the test, which is shown in Fig. 30, is very similar to that predicted with the semi-empirical 
method in the spreadsheet of Fig. 27. Figure 3 1 shows the notch in the acceleration input which 
resulted from the force limiting in Fig. 30. Comparison of force limited and unlimited data from 
low level runs indicates that the notch in the acceleration spectrum which results from force limiting 
is typically the mirror image of that portion of the force spectrum which would have exceeded the 
force limit, if no limiting were implemented. Notice the asymmetric shape of the notch in Fig. 31. 

It is characteristic that the notch, as well as the unlimited force, are steeper functions of frequency 
above a resonance than below it. This is in contrast to manual notches which are usually designed 
to be symmetric. 

The ACE spacecraft instrumentation includes a Spacecraft Launch Acceleration Measurement 
(SLAM) data acquisition system to measure, record, and transmit dynamic data at launch of the 
spacecraft on a Delta launch vehicle. The SLAM instrumentation includes a channel for the high 
frequency (20 to 2000 Hz) acceleration measured normal to the spacecraft honeycomb panel near 
one of the twelve mounting feet of the CRIS instrument and also a channel for the total normal 
force measured under the twelve mounting feet of the CRIS instrument. (The CRIS instrument is 
mounted on twelve uni-axial force transducers, and the output of these transducers is summed.) 
Figs. 32 and 33 show the spectral densities of these acceleration and force data measured during 
the prototype level acoustic test of the ACE spacecraft at NASA Goddard Space Flight Center 
(GSFC). The CRIS acceleration spectrum in Fig. 32 measured in the spacecraft acoustic test is 
approximately an order of magnitude less than the CRIS random vibration test acceleration input 
spectrum in Fig. 3 1 . This is unusually conservative, particularly when it is considered that the 
acoustic test spectrum is probably a conservative envelope of the flight acoustic levels. The CRIS 
normal force spectrum in Fig. 33 measured in the spacecraft acoustic test is approximately two 
orders of magnitude less than the CRIS random vibration test force limit spectrum in Fig. 30. One 
order of magnitude of conservatism in the force specification can be attributed to the conservatism 
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of the acceleration spectrum, since the force specification is proportional to the acceleration 
specification. However, the other order of magnitude of conservatism in the force spectrum must 
be attributed to the methods used to derive the force spectrum. It must be concluded that, even with 
force limiting, the random vibration test of the CRIS instrument was a severe overtest. 
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FIGURE 25. Ratio of Lateral External Force To Lateral Acceleration at Nominal CG for Cassini 
RPWS Instrument in DTM Spacecraft Acoustic Test 
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FIGURE 26. CRIS Instrument from Advanced Composition Explorer (ACE) Spacecraft 
Mounted with Flight Force Gages on Shaker for Vertical Vibration Test 
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FIGURE 27. Example of Spread Sheet for Calculating Force Limits- ACE Spacecraft CRIS 
Instrument 
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FIGURE 28. Tap Test Measurement of the Apparent Mass of the ACE Spacecraft Honeycomb 
Panel at One of the CRIS Instrument Mounting Positions 


Y Force 



FIGURE 29. Shaker Measurement of the Apparent Mass of the CRIS Instrument in 0.25 G Sine 
Sweep Test ( Multiply ordinate by four to obtain apparent mass.) 
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FIGURE 30. Total Vertical Force Limited to Specification in Protoflight Level Vertical Random 
Vibration Test of CRIS Instrument 



FIGURE 3 1 . Notched Acceleration Input in Protoflight Level Vertical Random Vibration Test of 
CRIS Instrument 
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FIGURE 32. Acceleration Spectrum Measured Normal to Panel Near Mounting Foot of CRIS 
Instrument in Protoflight Level Acoustic Test of ACE Spacecraft 
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FIGURE 33. Total Normal Force Spectrum Measured Under Twelve Mounting Feet of CRIS 
Instrument in Protoflight Level Acoustic Test of ACE Spacecraft 
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5.0 Instrumentation And Testing 


Herein are described the characteristics and use of piezo-electric force transducers and other 
instrumentation employed in force limited vibration testing. Some important considerations in the 
planning and conduct of the tests are also discussed. 


5.1 Piezo-electric Force Transducers 

The use of piezo-electric force transducers for force limited vibration testing is highly 
recommended over other types of force measurement means such as strain transducers, armature 
current, weighted accelerometers, etc. The historical review in Section 2 teaches that the basic 
concepts and methods of force limited vibration testing were espoused and recommended, some 
twenty or thirty years ago. It is the authors belief that the primary reason that force limiting has not 
been previously accepted and implemented on a routine basis was the lack of a practical force 
measurement device, prior to the development of the piezo-electric, tri-axial force transducer. The 
advent of these transducers has made the measurement of force in vibration tests almost as 
convenient and accurate as the measurement of acceleration. 

The high degree of linearity, dynamic range, rigidity, and stability of quartz make it an excellent 
piezo-electric material for both accelerometers and force transducers [16]. Similar signal 
processing, charge amplifiers and voltage amplifiers, may be used for piezo-electric force 
transducers and accelerometers. However, there are several important differences between these 
two types of measurement. Force transducers must be inserted between (in series with) the test 
item and shaker and therefore they require special fixtures, whereas accelerometers are placed upon 
( in parallel with) the test item or shaker. The total force into the test item from several transducers 
placed at each shaker attachment may be obtained by simply using a junction to add the charges 
before they are converted to voltage. On the other hand, the output of several accelerometers is 
typically averaged rather than summed. Finally, piezo-electric force transducers tend to put out 
more charge than piezo-electric accelerometers because the force transducer crystals experience 
higher loading forces, so sometimes it is necessary to use a charge attenuator between the force 
transducer and the charge amplifier. 


5.1.1 Force Transducer Preload 

Piezo-electric force transducers must be preloaded so that the transducer always operates in 
compression. The transverse forces are carried through the force transducer by friction forces. 
These transverse forces act internally between the quartz disks inside the transducer as well as 
between the exterior steel disks and the mating surfaces. Typically the maximum transverse load is 
0.2, the coefficient of friction, times the compressive preload. Having a high preload, and smooth 
transducer and mating surfaces, also minimizes several common types of transducer measurement 
errors, e.g. bending moments being falsely sensed as tension/compression if gapping occurs at the 
edges of the transducer faces. However, using flight hardware and fasteners, it is usually 
impossible to achieve the manufacturers recommended preload, so some calculations are necessary 
to insure proper performance. Sometimes it is necessary to trade-off transducer capability for 
preload and dynamic load. (This is often the case if there are large dynamic moments which can t 
be eliminated by designing the fixtures to align the load paths.) The three requirements for selecting 
the preload are: 1. it must be sufficient to carry the transverse loads through the transducer by 
friction, 2. it must be sufficient to prevent loss of compressive preload at any point on the 
transducer faces due to the dynamic forces and moments, and 3. it must be limited so that the 
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maximum stress on the transducer does not exceed that associated with the manufacturer’s 
recommended maximum load configuration. 


Transducer preloading is applied using a threaded bolt or stud which passes through the inside 
diameter of the transducer. With this installation, the bolt or stud acts to shunt past the transducer a 
small portion of any subsequently applied load, thereby effectively reducing the transducer’s 
sensitivity. Calibration data for the installed transducers is available from the manufacturer if they 
are installed with the manufacturer’s standard mounting hardware. Otherwise, the transducers must 
be calibrated in situ as discussed in the next section. 


5.1.2 Force Transducer Calibration 

The force transducer manufacturer provides a nominal calibration for each transducer, but the 
sensitivity of installed units depends on the size and installation of the bolt used for preloading and 
therefore must be calculated or measured in situ. This may be accomplished either quasi-statically 
or dynamically. Using the transducer manufacturer’s charge amplifiers and a low noise cable, the 
transducers will hold their charge for many hours, so that it is possible to calibrate them statically 
with weights or with a hydraulic loading machine. If weights are used, it is recommended that the 
calibration be performed by loading the transducers, re-setting to short-out the charge, and then 
removing the load, in order to minimize the transient overshoot. 

The simplest method of calibrating the transducers for a force limited vibration test is to conduct a 
preliminary low-level sine sweep or random run and to compare the apparent mass measured at 
low frequencies with the total mass of the test item. The appropriate apparent mass is the ratio of 
total force in the shaker direction to the input acceleration. The comparison must be made at 
frequencies much lower than the first resonance frequency of the test item. Typically the measured 
force will be approximately 80 to 90% of the weight in the axial direction and 90 to 95% of the 
weight in the lateral directions, where the preloading bolts are in bending rather than in tension or 
compression. Alternately, the calibration correction factor due to the transducer preloading bolt 
load path may be calculated by partitioning the load through the two parallel load paths according to 
their stiffness; the transducer stiffness is provided by the manufacturer, and the preload bolt 
stiffness in tension and compression or bending must be calculated. (The compliance of any 
structure in the load path between the bolt and transducer must be added to the transducer 
compliance.) 


5.1.3 Force Transducer Signal Conditioning 

It is strongly recommended that the total force in the shaker excitation direction be measured in a 
force limited vibration test. The total force from a number of transducers in parallel is readily 
obtained using a junction box which effectively sums the charges, and therefore the forces, before 
conditioning the signal with a charge amplifier. (An alternative is to specify limits for the force at 
individual attachment positions as in the case history in Section 6.1.2.) The same charge amplifiers 
used for piezo-electric accelerometers may be used for force transducers. However, the charge 
amplifiers made expressly for force transducers offer a choice of time constants, so that quasi-static 
(time constants of many hours) measurements of bolt preload, etc. may be obtained, as well as the 
dynamic measurements one usually associates with piezo-electric transducers. Also the charge 
amplifiers made expressly for force transducers usually have the capability to accommodate higher 
charge inputs, which are characteristic of force transducers. However, charge attenuators are 
readily available if they are needed. 
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Since vibration tests are normally conducted sequentially in three perpendicular axes, it is 
convenient to employ tri-axial force transducers. In addition, it is sometimes necessary to limit 
cross-axis force and the moments in addition to the in-axis force; this is particularly the case in tests 
of large eccentric test items such as spacecraft. For these applications, the six force msultant forces 
and moments for a single node may be measured with a combination commonly four, of tn-axia 1 
force transducers and commercially available voltage summers manufactured expressly for this 

purpose. 


5.2 Test Fixtures 


5.2.1 Design Concepts 

The preferred method of configuring the force transducers is to sandwich one transducer between 
the test item and conventional test fixture at each attachment position and use fasteners which are 
longer than the conventional ones to accommodate the height of the transducers. In this 
“S“„„ there is no fixture weight above the transducers and the transducer force ts identtcaj 

to the 8 force into the test item. Sometimes the preferred approach is jf 
many attachment points or the attachments involve shear pins m addition to bolts. In these cases 
mavbe necessary to use one or more light-weight intermediate adapter plates as an interface 
SSn the Sem and the force transducers. For example, if the test item mounts at three feet 
and each foot involves two bolts and a shear pin, a candidate design would be to have a small pla 
attached to a big stud for each foot. The small plate would pick up the two mounting bolts and 
shear pin, and the stud would go through a medium sized force transducer into a shaker ada P 
plate. Alternately, if the mounting configuration involves sixteen small bolts in a circular p , 
the fixture might consist of one intermediate ring which accepts the sixteen small bolts and 
mounted on eight equally spaced force transducers 

5.2.2 Fixture Weight Guideline 

The recommendation is that the total weight of any intermediate adapter plates above the force 
transducers do not exceed ten percent of the weight of the test item. This limitation is necessary 
because the force transducers read the sum of the force required to accelerate the interface plate and 
that delivered to the test item. If the fixture weight exceeds the 10% criterion, force limiting wi 
only be useful for the first one or two modes in each axis. Use of a circuit to subtract the mterfa 
plate force in real time, is not recommended because of the errors that result when the interface 
plate is not rigid. The use of armature current to measure shaker force 

because the weight of the armature and fixture typically are much greater than 10% of that ot th 
test item. 


5.2.3 Mass Cancellation 

If an intermediate plate is used and the plate moves unilaterally as a rigid body its acceleration may 
be measured with an accelerometer and subtracted, in real time from the gross f ° rce ^ e ^^ y 
the transducers underneath the intermediate plate to obtain the net force delivered to the test item. 
This approach has been utilized in a number of relatively simple tests, but is not recommended 
becauseof the possibility and consequences of errors. As the frequency approaches a resonance 
frequency of the intermediate plate on its mounting, the phase angle of the plate acceleration and 
applied force changes such that the aforementioned cancellation scheme does not work and may 
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make matters worse. Also, the additional instrumentation needed to implement the mass 
cancellation scheme gives rise to increased electrical noise and possibility of set-up error 


5.3 Testing Considerations 


5.3.1 Criteria For Force Limiting 

The purpose of force limiting is to reduce the response of the test item at its resonances on the 
shaker in order to replicate the response at the combined system resonances in the flight mounting 
configuration. Force limiting is most useful for structure-like test items which exhibit distinct, 
lg tly damped resonances on the shaker. Examples are complete spacecraft, cantilevered structures 
like telescopes and antennas, lightly damped assemblies such as cold stages, fragile optical 
components, and equipment with pronounced fundamental modes such as a rigid structure with 
flexible feet. The amount of relief available from force limiting is greatest when the structural 
impedance (effective mass) of the test item is equal to, or greater than, that of the mounting 
structure. However, it is recommended that notches deeper than 14 dB not be implemented without 
appropriate peer review. Force limiting is most beneficial when the penalties of an artificial test 
failure are high. Sometimes this is after an initial test failure in a screening type of test. 


5.3.2 Test Planning 

Several considerations need to be addressed in the test planning if force limiting is to be employed. 
First the size, number, and availability of the force transducers need to be identified as well as any 
special fixture requirements to accommodate the transducers. Next, the approach for deriving and 
updating the force specification needs to be decided. Finally the control strategy must be decided 
and written into the test plan. Special cases may include cross-axis force, moment, individual 
force, and response limiting in addition to or in lieu of the in-axis force. In some instances, the 
control strategy will be limited by the control system capabilities. In all cases, it is recommended 
that the control strategy be kept as simple as possible, in order to expedite the test and to minimize 
the possibility of mistakes. 

Accelerometers on the fixture are also required in force limited vibration tests in order to control the 
acceleration input to the acceleration specification at frequencies other than at the test item 
resonances. In addition, it is often convenient to use a limited number of accelerometers to measure 
the response at critical positions on the test item. These response accelerometers may be used only 

for monitoring or, if justified by appropriate rationale, for response limiting in addition to the force 
limiting. 


5.3.3 Cost And Schedule 


Project managers often inquire regarding the cost and schedule impact of doing force limiting. 
Once force limiting has become routine in the vibration testing laboratory, its use definitely results 
in both cost and schedule savings. The primary savings is through the prevention of unnecessary, 
overtest failures, which adversely impact both cost and schedule. Implementation of force limiting 
on a routine basis also eliminates most of the contentious discussions about test levels. Also the 
effort previously required to measure, analyze, and limit numerous responses in complex vibration 
tests is greatly reduced by force limiting, as illustrated by the Cassini spacecraft test example in 
Section 6.2. v 
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The first or second time that an organization employs force limiting, there will naturally be some 
down time and slow going. The first consideration must be the availability ot force transducers of 
the appropriate size. Tri-axial force transducers are relatively expensive, compared to 
accelerometers, and are sometimes long-lead procurement items, so an assortment of transducers 
are usually procured over time by several projects and maintained by the vibration test laboratory 
for use on future projects. If the first application is a flight project, it is recommended that some 
help be sought from the sponsor or another organization that has experience in force limiting. The 
only other potential additional cost of using force limiting is the development of special fixtures 
Usually the force transducers can be utilized simply by placing one at each mounting position and 
using a longer bolt to accommodate the thickness of the transducer. Configurations where special 
fixtures may be needed are: those which involve a large number of mounting points, say more than 
twelve, or those with shear pins or complicated fittings, such as the flight latches in the example of 
Section 6.1.2. 


5.3.4 Specification of Force Limits 

Force limits are analogous and complementary to the acceleration specifications used in 
conventional vibration testing. Just as the acceleration specification is the frequency spectrum 
envelope of the in-flight acceleration at the interface between the test item and flight mounting 
structure, the force limit is the envelope of the in-flight force at the interface. In force limited 
vibration tests, both the acceleration and force specifications are needed, and the force specification 
is proportional to the acceleration specification. Therefore force limiting does not compensate for 
errors in the development of the acceleration specification , e.g. undue conservatism or lack thereof. 
These errors will carry over into the force specification. Since in-flight vibratory force data are 
lacking, force limits are usually derived from coupled system analyses and impedance information 
obtained from measurements or finite element models (FEM). Also, considerable data on the 
interface force between spacecraft and components are becoming available from spacecraft acoustic 
tests, and semi-empirical methods of predicting force limits are available. 

Force spectra have typically been developed in one-third octave bands (see example in Section 
4.4), but other bandwidths, e.g. octave or one-tenth octave bands, may also be used. Force 
limiting may usually be restricted to an upper frequency encompassing approximately the first three 
modes in each axis; which might be approximately 100 Hz for a large spacecraft, 500 Hz for an 
instrument, or 2000 Hz for a small component. It is important to take into account that the test item 
resonances on the shaker occur at considerably higher frequencies than in flight. Therefore care 
must be taken not to roll off the force specification at a frequency lower than the fundamental 
resonance on the shaker and not to roll off the specification too steeply, i.e. it is recommended that 
the roll-offs of the force spectrum be limited to approximately 9 dB/octave. 


5.3.5 Vibration Controllers 

Most of the current generation of vibration test controllers have the two capabilities needed to 
implement force limiting. First, the controller must be capable of extremal control, sometimes 
called maximum or peak control by different vendors. In extremal control, the largest of a set of 
signals is limited to the reference spectrum. (This is in contrast to the average control mode in 
which the average of a set of signals is compared to the reference signal.) Most controllers used in 
aerospace testing laboratories support the extremal control mode. The second capability required is 
that the controller must support different reference spectra for the response limiting channels, so 
that the force signals may have limit criteria specified as a function of frequency. Controllers which 
support different reference spectra for limit channels are now available from most venders and in 
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addition upgrade packages are available to retrofit some of the older controllers for this capability. 
If the controller does not have these capabilities, notching of the acceleration specification to limit 
the measured force to the force specification must be done manually in low level runs. 


5.3.6 Low-Level Runs 

It is advantageous to keep the number of test runs as low as feasible, both to save testing time and 
to avoid accumulating unnecessary fatigue of the test hardware. A low-level sine-sweep or random 
run with a flat frequency spectrum is often conducted before and after the high-level vibration run 
in each axis to measure the vibration signature for “health” monitoring. The reaction force is an 
excellent choice for the health monitoring signature, and the ratio of force to input acceleration 
measured in the before” run can also be used to update the effective masses and resonance 
frequencies used to derive the force specification. (Sometimes the derivation of the force 
specification is deferred until this data become available.) 

It is also advantageous to conduct two low-level runs (often -18 dB) with the same input 
acceleration spectral shape as the high-level run; the first without force limiting and the second with 
force limiting, using a scaled down force limit. Comparison of the forces measured in these two 
low-level runs with the amount of notching achieved in the second run with force limiting provides 
verification that the force limit is appropriate and that the notching is as it should be. Extrapolation 
ot measured responses to full level and comparison with the loads criteria and any adjustments to 
the force specification should take place after these two runs. After it is determined that the results 
are satisfactory, any intermediate-level runs and the full-level run may be conducted. With the new 
controllers, it is becoming common practice to come “on the air” at a low level and then proceed 
through the intermediate levels to full level, without shutting down. Thus ideally, each axis may be 
conducted with no more than five runs: 1) low-level pre-test signature, 2) -18 dB without force 
limiting, 3) -18 dB with force limiting, 4) “on the air” at intermediate level and progressing to full 
level, and 5) low-level post-test signature. 



6.0 Case Histories 

The three random vibration tests selected as case histories include a component, an instrument, and 
a spacecraft. These test items span a mass range greater than 101 The component test was 
conducted in March of 1991 and the spacecraft test in November of 1996; the technology, 
particularly as regards the prediction of force limits and the shaker control, progressed 
considerably in this time interval. 


6.1 Hubble Wide Field Planetary Camera II AFM Component and 
Complete Instrument 


6.1.1 Articulating Fold Mirror Component Random Vibration Test [42] 

There are three articulating fold mirror (AFM) assemblies in the Wide Field Planetary Camera II 
(WFPCII) installed in the Hubble Telescope during the first servicing mission in December 1993. 
The role of the AFM is to provide a means for very accurate on-orbit alignment of the optical beam 
on the secondary relay mirrors which contain the correction for the Hubble primary mirror 
spherical aberration. The photograph in Fig. 34 shows an AFM before the mirror is coated. The 
AFM utilizes three small electro-strictive actuators to articulate the mirror. By necessity, the AFM 
is small and delicate; the total unit weighs approximately 100 gm and the articulating portion, the 
mirror and bezel, weighs less than 30 gm. 

Figure 35 shows an AFM mounted in the vibration test fixture designed to accommodate three 
small (approx. 2.5 cm square) commercially available tri-axial force transducers. Conical spacers 
are used to mate the transducers with the #2 screws which normally attach the AFM to a bulkhead 
of the optical bench of the WFPCII. The apparent mass normal to the WFPCII bulkhead was 
measured with an impact hammer to be 3.2 kg at 800 Hz, the first axial resonance of the AFM. The 
force specification was determined using the simple TDFS method. Notice from Fig. 5 that even 
with the small ratio (0.01) of AFM mass to bulkhead effective mass and with the measured Q of 
50, the normalized ratio of force spectral density to acceleration spectral density of (100) is still 
14 dB less than Q squared (2500), which is the amplification expected in a conventional vibration 
test. 

The results of the Z-axis protoflight random vibration test of the AFM qualification unit are shown 
in Fig. 36. The unnotched input acceleration spectral value of 0.004 GVHz at 800 Hz was 
determined from measurements during the acoustic test of the optical bench of WFPCI, the original 
camera on the Hubble Space Telescope. Figure 36 also shows the total axial force input to the 
AFM, i.e. the sum of the Z axis outputs of the three force transducers in Fig. 35. The force shown 
in Fig. 36 is actually scaled up from a low level (18 dB down) test without force limiting to show 
what the force would have been in the full level test if the force were not limited. In the signal 
conditioning, the force feedback signal was multiplied by the ratio of the acceleration specification 
to the force specification so that extremal control of the both force and acceleration could be 
implemented by comparing both signals to the acceleration specification, as discussed in Section 
5.3.2. For this reason, the force signal in Fig. 36 can be compared directly with the acceleration 
specification. Figure 36 shows that in the protoflight-level test with force limiting, the controller 
automatically notched the acceleration input by the amount the unlimited force signal would have 
exceeded its specification, i.e. about 10 dB, which is 4 dB less than that estimated from Fig. 5. 
The notch is very sharp and approximately the mirror image of the force peak. It is impractical to 
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manually put in such sharp notches. Also without the force sensors to detect the frequency of the 
force peak, it would be difficult to place the notch at the correct frequency. 


6.1.2 Wide-Field Planetary Camera II Instrument Random Vibration Test [42] 


The complete WFPCII (weight 284 kg) was subjected to a vertical axis protoflight level random 
vibration test as shown in Fig. 37. A large (approx. 10 cm diam.) commercially available, tri-axial 
force transducer was located just below each of the three latches, which fasten the camera to the 
Hubble telescope. Force specifications for all three directions at all three latches were derived using 
the simple TDFS method with apparent mass data for the WFPCII and for the honeycomb 
container used to transport the WFPCII in the Space Shuttle to rendezvous with the orbiting 
telescope. The apparent masses of the container were measured by NASA GSFC with a modal 
impact hammer. The data for the vertical direction at the A latch of the container is shown in Fig. 
38. (The A latch is on the far left in Fig. 37.) The three components (X, Y & Z) of force at each of 
the three latches were recorded and analyzed in a low level sine test of the WFPCII preceding the 
random vibration test, and the apparent mass of the WFPCII in the vertical direction at the A latch 
is shown in Fig. 39. After a low-level random run, it was determined that three of the four control 
channels available on the older controller were essential to control the high frequency acceleration 
input at each of the three latches, so only one control channel was available for force limiting. 

(New controllers have 32 to 64 control channels.) On the basis of analysis and the low level data, it 
was decided to limit the vertical force at the A latch which reacts most of the load because of the 
outboard (to the right in Fig. 37) center-of-gravity of the camera. 

The acceleration measured at each of the three latches is compared with the acceleration 
specification for the protoflight random vibration test in Fig. 40, and the vertical force measured at 
each of the three latches is compared to the force specification for the A latch in Fig. 41. 
Comparison of Figs. 40 and 41 shows that limiting the A latch vertical force resulted in notching of 
the acceleration input at 35 Hz, 100 Hz and 150 Hz. The acceleration notch at 70 Hz was effected 
by reducing the vertical force specification at 70 Hz to compensate for the A latch transverse force 
not being in the control loop. Above 300 Hz, one of the three control accelerometers is equal to the 
specification at every frequency. Comparison of the notched vibration input in the AFM component 
vibration test described in the last section with the corresponding vibration levels measured a year 
later on the optical bench in the WFPCII system acoustic test confirmed that the vibration levels in 
the component test enveloped those in the system acoustic test, as they should for a valid 
component screening test. 


6.2 Cassini Spacecraft System Random Vibration Test [28, 43] 

Figure 42 is a photograph of the Cassini flight spacecraft mounted on the shaker for the vertical 
random vibration test which was conducted in November of 1997 [28, 43]. The Cassini orbiter 
weighs a total of 2,150 kilograms (4,750 pounds); after attaching the 350-kilogram Huygens probe 
(on the right in Fig. 42) and a launch vehicle adapter and loading more than 3,000 kilograms 
(6,600 pounds) of propellants, the spacecraft weight at launch is about 5,800 kilograms (12,800 
pounds). The weight for the vibration test was somewhat less 3,809 kg (8,380 lb.), because the 
tanks were loaded to 60% of capacity with referee fluids. Because of the very dim sunlight at 
Saturn’s orbit, solar arrays are not feasible and power will be supplied by a set of three 
Radioisotope Thermoelectric Generators (RTG’s) which use heat from the natural decay of Most of 
plutonium to generate electricity to run Cassini. (Two of the RTG’s are visible in the lower center 
and left of Fig. 42.) Twelve science experiments are carried onboard the Cassini orbiter and 
another six fly on the Huygens Titan probe. The schematic in Fig. 7 indicates a number of the 
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Cassini spacecraft instruments, the orbiter instruments are mounted on the Remote Sensing 
Platform at the upper left and on the Fields and Particles Platform at the upper right in Fig. 42. 

Figure 43 shows the plan view of the spacecraft mounting ring before the spacecraft is attached to 
the shaker The black offset weight positioned in the upper right quadrant of the ring is being 
subjected to vibration for moment proof testing of the shaker and mounting configuration. Ihe 
spacecraft bolts to the ring at eight positions corresponding to the mounting feet locations on the 
spacecraft/ launch vehicle adapter. A large tri-axial force transducer is located under the mounting 
ring at each of these eight positions. The shaker fixture is restrained from moving laterally during 
the & spacecraft vertical vibration test by three hydraulic bearings. The force capability of the shaker 
is about 35,000 lb. and virtually all of this capability was used to vibrate the spacecraft and shaker 
fixtures, which weighed about 6,000 lb. in addition to the 8380 lb. spacecraft. 


Figure 44 compares the vibration test acceleration input specification with launch vehicle 
specifications and with data from a previous Titan launch vehicle flight. TJie acceleration 
specification was originally somewhat higher (0.04 G /Hz compared to 0.01 G/Hz). The 
specification was lowered in the 10 to 100 Hz frequency regime after reviewing the results of an 
extensive FEM pre-test analysis, which indicated that excessive notching would be required with 
the higher-level input. The specification was subsequently lowered in the 100 to 200 Hz regime as 
well in order to accommodate the force capability of the shaker power amplifier, which was over 
ten years old and exhibited some instability problems during the two month period preceding the 
Cassini spacecraft vibration test. The resulting 0.01 GVHz specification is less than the Booster 
Powered Phase specification at frequencies greater than approximately 80 Hz, but exceeds the 
Maximum Envelope of the TIV-07 Flight Data. The acceleration specification is defined in 

Fig. 45. 


Figure 46 shows the force specification for the Cassini flight spacecraft random vibration test. 1 he 
specification was derived by multiplying the acceleration specification in Fig. 45 by the squared 
weight of the spacecraft and by a factor of one-half. (This corresponds to the semi-empirical 
method discussed in Section 4.2. 1 with a value of C = 0.707.) This value of C was selected on the 
basis of the pre-test analysis and in order to keep the proof test, which had a margin of 1.25 over 
the test limit loads, within the shaker force capability. The force specification was not rolled off at 
the shaker fundamental resonance as shown in Equation 29, because neither the pre-test analysis 
nor the actual vibration test data exhibited a distinct fundamental resonance of the spacecraft in the 
vertical axis. During the test, it was not necessary to modify or update the force specification 
specified in the test procedure. 


Figures 47 and 48 respectively, show the input acceleration and force spectra measured in the 
actual full-level vibration test. Comparison of the measured acceleration spectra with the 
specification in Fig. 47 shows significant notching of ~8 dB at the probe resonance of 
approximately 17 Hz and of ~14 dB at the tank resonance of approximately 38 Hz. The responses 
at a number of critical positions on the spacecraft, as well as the other five components of the total 
input force vector, were monitored during the testing, but only the total vertical force signal was 
used in the controller feedback to notch the acceleration input. Comparison of the measured force 
with the specified force in Fig. 48 verifies that the force was at its limit at all the frequencies where 
notching occurred in the input acceleration. 


Figures 49 and 50 show the acceleration inputs measured near the feet of a number of instruments 
mounted on the Fields and Particles and Remote Sensing Pallets, respectively. Comparison of 
these measured data with the random vibration test specifications for the instruments, which are 
also indicated in Figs. 49 and 50, demonstrates that many of the instruments reached their 
component vibration test levels in the spacecraft vibration test. The significant excedances below 
50 Hz are covered by the instrument sine vibration test equipment. In addition, several major 
components of the spacecraft including the Huygens probe upper strut, the three RTG s, the 
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magnetic canister struts, and the Fields and Particles Pallet struts reached their flight limit loads 
during the spacecraft vibration test. The only anomaly after the test, other than possibly those ' 
associated with spacecraft functional tests for which data are not available, was that the electrical 
resistance between the engineering model RTG and the spacecraft structure was measured after the 
test and found to be less than specified. The insulation between the RTG adapter bracket and the 
spacecraft was redesigned to correct this problem. 



FIGURE 34. Photograph of 30-Gram Articulating Fold Mirror (AFM) Before Coating 



Y-Axis 


FIGURE 35. Sketch of Vibration Test Fixture for AFM Showing Mounting on Three Small 
Tri-axial Force Gages Using Conical Adapters 
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FIGURE 36. Notched Acceleration Input in Vertical Protoflight Random Vibration Test of the 
WFPCI1 Articulating Fold Mirror (AFM) 



FIGURE 37. Hubble Space Telescope 600-Kilogram Wide Field Planatary Camera II (WFPCII) 
Mounted on Shaker for Vertical Axis Vibration Test 
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Dynamic Mass (kg.). 



FIGURE 39. Apparant Mass of WFPCII at A Latch in Vertical Direction Measured on Shaker in 
Low-Level Sine-Sweep Vibration Test 
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FIGURE 40. Measured and Specified Acceleration Input Spectra at Three Latches of the Wide 
Field Planetary Camera II in the Protoflight Random Vibration Test 



10 100 1000 
Frequency, Hz 

FIGURE 41. Measured and Specified Axial Force Spectra at Three Latches of the Wide Field 
Planetary Camera II in the Protoflight Random Vibration Test 
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FIGURE 42. Cassini Spacecraft Mounted on Shaker for Vertical Random Vibration Test 
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FIGURE 43. Plan View of Spacecraft Mounting Ring with Offset Weight for Moment Proof Test 
and of Lateral Restraint System for Cassini Spacecraft Vibration Test 


6-10 


u 


(2H ciSd 



FIGURE 44. Comparison of Cassini Spacecraft Random Vibration Acceleration Input 
Specification with Launch Vehicle Specifications and Flight Data 
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FIGURE 46. Cassini Flight Spacecraft Random Vibration Test Force Specification 



FIGURE 47 Cassini Flight Spacecraft Random Vibration Test Full-Level Input Acceleration 
Spectral Density (Notches from Specification Due to Force Limiting) 
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FIGURE 48. Total Vertical Force Measured in Cassini Flight Spacecraft Full-Level Random 
Vibration Test (Comparison with Vertical Force Limit Specification) 
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FIGURE 49. Acceleration Inputs to Fields and Particles Pallet Instruments in Cassini Full-Level 
Random Vibration Test (Comparison with Instrument Random Vibration Test 
Specification) 
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FIGURE 50. Acceleration Inputs to Remote Sensing Pallet Instruments in Cassini Full-Level 
Random Vibration Test (Comparison with Instrument Random Vibration Test 
Specification) 
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7.0 Conclusions 


Here the key points discussed in this monograph .ire reiterated, and some suggestions for 
further development of the technology are offered. 


1. Technical literature reviewed supports the conclusion that the pioneers of aerospace 
vibration testing recognized the dangers of vibration overtesting and understood that the 
very high shaker impedance was the culprit f 1.2]. Many researchers during the 
subsequent thirty years studied the mounting structure impedance problem and 
developed conceptual solutions, many of which are a part of the force limiting approach 
described herein [4,10], One can only conclude that either the instrumentation or the 
need was not sufficient to drive the implementation of this technology until recently. It 
is the author s belief that the recent advent of piezo-electric, tri-axiai force transducers 
was the enabling factor. 

2. In built-up aerospace configurations, the structural impedance of the mounting structure 
and test item are comparable and there is little amplification at the resonance frequencies 
of the test item. The interface acceleration has notches at the test item resonance 
frequencies, due to the vibration absorber effect. 


3. The shaker, by contrast, has very high mechanical impedance and the test item can have 
very large amplification at its resonance frequencies, in addition, the test item resonance 
frequencies on the snaker occur at significantly higher frequencies than the coupled 
ystem resonances in the built-up configuration. 


4. The object of force limited vibration testing is to make the input force at the test item 
resonances in the vibration test equal to the maximum interface force in the flight 
configuration. The goal is to replicate the internal forces and stresses in the flight 
environment, but the simulation is not exact because the resonance frequencies and 
mode shapes are different. The result of limiting the input force at resonances in the 
vibration test is that the input acceleration is notched in a manner similar to that due to 
the vibration absorber effect in flight. However, the notch frequencies will be slightly 
different on the shaker and in flight, due to the differences in off-axis boundary 
conditions. 


5. The effective mass concept, developed in the early 70’s 115]. provides a theoretical 
basis for analyzing structural impedance data and FEM results. Until recently, the 
effective mass concept was know and used by only a select group of analysts, but the 
concept is gaming acceptance and it is hoped that an increased acceptance of the concept 
will be facilitated by the discussion in Section 3.3.4. provided by one of the concept's 
originators, as to the derivation of the effective mass from NASTRAN of other FEMs. 


6. The philosophy adopted in the derivation of force specifications is to start with the 
traditional acceleration specification, which is the envelope of the flight interface 
acceleration, and add an analogous force specification, which is the envelope of the 
flight interface forces. Thus the force specifications derived analytically herein are 
proportional to the assumed acceleration specification, and any errors or undue 
conservatism in the acceleration specification carry over into the force specification. 
Thus force limiting should not be perceived as a method of compensating for errors in 
the acceleration specification. Rather, it is a method of automatically inserting notches 
in the acceleration spectrum at the proper frequencies and of the proper depth. 



7. The general method of deriving force specifications is to develop a coupled system 
model of the source and load, with the modal parameters of each determined from FEM 
analysis or from impedance measurements. Then the ratio of the frequency envelope of 
the interface force to the envelop e of the interface acceleration is determined. Finally, 
this ratio is multiplied by the acceleration specification, developed in a conventional 
manner (usually semi-empirically), to obtain the force specification. Two applications 
of this general method involving a simple and a complex TDFS are derived herein and 
the results are presented parametrically [27]. Other applications of the method are 
available in the literature [24] and it is envisioned that more sophisticated models will be 
developed in the future. 

8. The literature [3,10,14] and spacecraft system acoustic tests [28] have provided data for 
the development of semi-empirical force limits which are much simpler to apply and 
appear to yield satisfactory results. Little flight data are currently available for vibratory 
force, but several flight measurement programs are in progress. 

9. The advent of tri-axial, piezo-electric force transducers and of a new generation of 
digital controllers have facilitated the application of force limited vibration testing. The 
piezo-electric force transducers are easy-to-use, rugged, compact, have a wide dynamic 
range, and can readily be configured to measure all six components of force and 
moment [16]. The convenient measurement of total external force in a vibration test 
now makes it possible to measure the acceleration of the center-of-gravity, so that the 
design capability of aerospace structures can be conveniently verified in a vibration test. 

10. There are many topics in the development of force limited vibration testing which 
require further investigation. Impedance methods are most convenient when two 
systems are connected at a single node, and this assumption is implicit herein. Many 
difficulties occur when one considers multi-point connections [8], which are almost 
always the case in the real world. The issue of overtesting due to uncorrelated inputs at 
multiple attachment points has not been addressed herein. Future vibration controllers 
will probably offer the capability to control phase, and then appropriate phase 
specifications between force and acceleration will have to be developed [4, 13]. New 
force transducers are being developed, including devices which can both measure and 
generate force, and new instrumentation developments will open the door to new 
testing techniques. Comparison of flight data with the prediction methods discussed 
herein will almost certainly give rise to some discrepancies, which will need to be 
resolved. 


7-2 


1 



Appendix A Definition of Symbols 


A = interface acceleration 

A b = base acceleration 

A 0 = free acceleration of source 

A s = acceleration specification 

C = dashpot constant 

C = constant 

F = interface force 

F u = blocked force of source 

F 5 = force specification or limit 

F e = excitation force 

k = spring stiffness 

k = physical stiffness matrix 

M 0 = total mass 

M = residual mass 

m = modal mass 

M = apparent mass, F/A 

m = physical mass matrix 

M = modal mass matrix 

Q = dynamic amplification factor 

S AA = acceleration spectral density 

Spp = force spectral density 

u = absolute displacements 

U = generalized modal displacement 

a = ratio of modal to residual masses 

P = ratio of analysis frequency to resonance frequency 

(]) = mode shape 

(i = ratio of load to source residual masses 
(D = radian frequency 

co o = natural frequency of uncoupled oscillator 

£2 = ratio of load to source uncoupled resonance frequencies 

Subscripts 

1 = source oscillator 

2 = load oscillator 

F = unrestrained (free) 

P = prescribed 

N = modal set 

R = rigid body set 

n = single mode 

p = reaction force direction 

q = prescribed acceleration direction 
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FORWARD 


This handbook is approved for use by NASA Headquarters and all Field centers and is intended to 
provide a common framework for consistent practices across NASA programs. 

The primary goal of vibration tests of aerospace hardware is to identify problems which, if not 
remedied, would result in flight failures. This goal can only be met by implementing a realistic (flight- 
like) test with a specified positive margin. In most cases, the goal is not well served by traditional 
acceleration controlled vibration tests which indeed screen out flight failures but in addition, cause 
failures, which would not occur in flight. This overtest or “screening” test approach, which may 
have served its purpose in the past, is too expensive and inefficient for today’s environment of low- 
cost missions. The penalty of overtesting is manifested in design and performance compromises, as 
well as in the high costs and schedule overruns associated with recovering from artificial test failures. 

It has been known for thirty years that the major cause of overtesting in aerospace vibration tests is 
associated with the infinite mechanical impedance of the shaker and the standard practice of 
controlling the input acceleration to the frequency envelope of the flight data. This approach results in 
artificially high shaker forces and responses at the resonance frequencies of the test item. To alleviate 
this problem it has become common practice to notch the input acceleration to limit the responses in 
the test to those predicted for flight, but this approach is very dependent on the analysis, which the 
test is supposed to validate. Another difficulty with response limiting is that it requires placing 
accelerometers on the test item at all the critical locations, many of which are often inaccessible, and 
which in the case of large test items, involves extensive instrumentation. 

The advent of new instrumentation has made possible an alternative, improved vibration testing 
approach based on measuring and limiting the reaction force between the shaker and test item. The 
major break through is the availability of piezo-electric triaxial force gages developed for other 
commercial markets. Piezo-electric force gages are robust, relatively easy to install between the test 
item and shaker, and require the same signal conditioning as piezo-electric accelerometers commonly 
used in vibration testing. Also, a new generation of vibration test controllers now provide the 
capability to limit the measured forces and thereby notch the input acceleration, in real time. To take 
advantage of this new capability to measure and control shaker forces, a rationale for predicting the 
flight limit forces has been developed and applied to many flight projects during the past five years. 
Force limited vibration tests are now conducted routinely at the Jet Propulsion Laboratory (JPL) and 
also at several other NASA centers, government laboratories, and many aerospace contractors. 

This handbook describes an approach which may be used to facilitate and maximize the benefits of 
applying this relatively new technology throughout NASA in a consistent manner. A monograph, 
which provides more detailed information on the same subject, is also available for reference. 

Requests for information, corrections, or additions to this handbook should be directed to the 
Mechanical Systems Engineering and Research Division, Section 352, Jet Propulsion Laboratory, 
4800 Oak Grove Dr., Pasadena, CA 91109. Requests for additional copies of this handbook should 
be sent to NASA Engineering Standards, EL02, MSFC, AL, 35812 (telephone 205-544-2448). 


Daniel R. Mulville 
Chief Engineer 
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1.1 Purpose 

This handbook establishes a methodology for conducting force limited vibration tests for all 
NASA flight projects. The purpose is to provide an approach which may be consistently 
followed by those desiring to use force limiting, without having to conduct an extensive 
literature search or research and development effort before conducting the test. The decision 
to use, or not use, force limiting on a specific project and in a specific vibration test, and 
the responsibility for applying the method correctly, are left to the project or the cognizant 
engineer. A monograph on Force Limited Vibration Testing is available for reference and it 
is recommended for those needing more detailed technical information. 

1 .2 A pplicability 

This handbook recommends engineering practices for NASA programs and projects. It 
may be cited in contracts and program documents as a technical requirement or as a 
reference for guidance. Determining the suitability of this handbook and its provisions is 
the responsibility of program/project management and the performing organization. 
Individual provisions of this handbook may be tailored (i.e., modified or deleted) by 
contract or program specifications to meet specific program/project needs and constraints. 

For the purpose of this handbook, a force limited vibration test is any vibration test in 
which the force between the test item and shaker is measured and controlled. (The 
recommended means of measuring the force is with piezo-electric force gages, but other 
means, e.g. shaker armature current or strain gages, may be useful in special situations. 
Similarly, the control of the force is preferably accomplished in real time, but iterative, off- 
line control may be employed as a stepping stone.) If the force is not measured and 
controlled, the test is not considered a force limited vibration test, and this handbook does 
not apply. This distinction is important because in the past some have found it convenient 
to simulate a force limited test and then to use the analytical results to notch the acceleration 
input in the test. The simulation approach is not recommended because measurement of the 
force is considered to be the essential element of the force limiting approach. 

The handbook is applicable to all force limited vibration tests of NASA flight hardware 
including launch vehicle and aircraft equipment, spacecraft, instruments, and components. 
However since the purpose of force limiting is to mitigate the effect of test item resonances 
in the vibration test, the technique is most useful for structure-like equipment and for fragile 
equipment such as optics and complex instruments. 

2.0 APPLICABLE DOCUMENTS 

2. 1 General . The applicable documents cited in this handbook are listed in this section only 
for reference. The specified technical requirements listed in the body of this document must 
be met whether or not the source document is listed in this section. 

2.2 Government documents . The following Government documents form a part of this 
document to the extent specified herein. Unless otherwise specified, the issuances in effect 
on date of invitation for bids or request for proposals shall apply. 

NASA - RP- 1403 - Force Limited Vibration Test 

Monograph, May 1997 

NASA - GUIDELINES - XXXX - Dynamic Environment 

Guidelines, September 1997 
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NASA -STD -7001 


Payload Vibroacoustic 
Test Criteria, June 21, 1996. 


NASA - STD - 7002 - Payload Test Requirements, 

July 10, 1996. 

2.3 Non-government publications . The following documents a part of this document to the 
extent specified herein. Unless otherwise specified, the issuances in effect on the date of 
invitation for bids or request for proposals shall apply. 

Blake R. E., “The Need to Control the Output Impedance of Vibration and Shock 
Machines”, Shock and Vibration and Associated Environments, Bulletin No. 23, 1954. 

Salter, J. P., “Taming the General-Purpose Vibration Test”, Shock and Vibration and 
Associated Environments, Bulletin No. 33, Part III, March 1964, pp. 211-217. 

Murfin, W. B., “Dual Specifications in Vibration Testing”, Shock and Vibration Bulletin, 
No. 38., Part 1, 1968, pp. 109-113. 

Wada, B. K., Bamford, R., and Garba, J. A., “Equivalent Spring-Mass System: A 

Physical Interpretation”, Shock and Vibration Bulletin, No. 42, 1972, pp. 215-225. 

Scharton, T. D., Boatman, D. J., and Kern, D. L., “Dual Control Vibration Testing”, 
Proceedings of 60th Shock and Vibration Symposium, Vol. IV, 1989, pp. 199-217. 

Smallwood, D. 0„ “An Analytical Study of a Vibration Test Method Using Extremal 
Control of Acceleration and Force”, Proceedings of Institute of Environmental Sciences 
35th Annual Technical Meeting, 1989, pp. 263-271. 

Scharton, T. D., “Analysis of Dual Control Vibration Testing”, Proceedings of Institute of 
Environmental Sciences 36th Annual Technical Meeting, 1990, pp. 140-146. 

Smallwood, D. O., “Development of the Force Envelope for an Acceleration/Force 
Extremal Controlled Vibration Test”, Proceedings of 61st Shock and Vibration 
Symposium, Vol. I, 1990, pp. 95-104. 

Scharton, T. D., “Force Limited Vibration Testing at JPL”, Proceedings of the Institute of 
Environmental Sciences 14th Aerospace Testing Seminar, 1993, pp. 241-251. 

Scharton, T. D., “Vibration-Test Force Limits Derived from Frequency-Shift Method”, 
AIAA Journal of Spacecraft and Rockets, Vol. 2, No. 2, March 1995, pp. 312-316. 

Scharton, T., Bamford, R., and Hendrickson, J., “Force Limiting Research and 
Development at JPL”, Spacecraft and Launch Vehicle Technical Information Meeting, 
Aerospace Corp., Los Angeles, CA, June 7, 1995. 

Chang, K. Y. and Scharton, T. D., “Verification of Force and Acceleration Specifications 
for Random Vibration Tests of Cassini Spacecraft Equipment”, ESA/CNES Conference on 
Spacecraft Structures, Materials, and Mechanical Testing, Noordwijk, NL, March 27-29, 
1996. 

2.4 Order of precedence. Where this document is adopted or imposed by contract on a 
program or project, the technical guidelines of this document take precedence, in the case of 
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conflict, over the technical guidelines cited in other referenced documents. This handbook 
does not apply to payload programs approved prior to the date of this document. Also, this 
handbook does not address safety considerations that are covered thoroughly in other 
documents; but if a conflict arises, safety shall always take precedence. Nothing in this 
document, however, supersedes applicable laws and regulations unless a specihc 
exemption has been obtained. 

3.0 DEFINITIONS 

ACCELERANCE— Complex frequency response function which is ratio of acceleration to 
force 

ACCELERATION OF C.G. -Acceleration of instantaneous centroid of distributed masses 
(equal to external force divided by total mass, according to Newton s 2nd Law) 

APPARENT MASS-Complex frequency response function which is ratio of force to 
acceleration 

CONTROL SYSTEM-The hardware and software which provides means for the test 
operator to translate vibration specifications into the drive signal for the shaker 

DESIGN VERIFICATION TEST-Test to see if as-built test-item can survive design loads 

DUAL CONTROL— Control of both force and vibration 

DYNAMIC ABSORBER— SDFS tuned to excitation frequency to provide reaction force 
which reduces motion at attachment point 

EFFECTIVE MASS— Masses in model consisting of SDFS’s connected in parallel to a 
common base, so as to represent the apparent mass of a base driven continuous system. 
The sum of the effective modal masses equals the total mass. 

EXTREMAL CONTROL— A shaker controller algorithm based on control of the maximum 
(extreme) of a number of inputs in each frequency control band 

FLIGHT LIMITS-Definition of accelerations or forces which are believed to be equal to 
the maximum flight environment, often P(95/50) 

FORCE LIMITING-Reduction of the reaction forces in a vibration test to specified values, 
usually to the interface forces predicted for flight, plus a desired margin 

IMPEDANCE-Complex frequency response function which is ratio of force to velocity 
quantities (Sometimes used to refer to ratio of force to any motion quantity.) 

LEVEL-Test input or response in decibels (dB), dB = 20 log amplitude = 10 log power 
LOAD- Vibration test item 

MARGIN-Factor to be multiplied times, or decibels to be added to, the flight limits to 
obtain the test specification 

NOTCHING-Reduction of acceleration input spectrum in narrow frequency bands, 
usually where test item has resonances 
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QUALITY FACTOR-Measure of the amplification of the response at resonance, equal to 
the reciprocal of twice the critical damping ratio. 1 

QUASI-STATIC ACCELERATION-Combination of static and low frequency loads into 
an equivalent static load specified for design purposes as C.G. acceleration 

RESIDUAL MASS— Sum of the effective masses of all the vibration modes with 
resonance frequencies greater than the excitation frequency. 

RESPONSE LIMITING— Reduction of input acceleration to maintain measured response 
at or below specified value, usually as predicted for flight plus desired margin 

SHAKER-The machine which provides vibratory motion to the test item, usually 
electrodynamic, in aerospace testing ( can also be hydraulic or rotary) 

THREE AXIS LOAD CELL-Force gage which measures the three perpendicular 
components of force simultaneously 

SINGLE DEGREE-OF-FREEDOM SYSTEM (SDFS) —Vibration model with one mass 

SOURCE-Test item support structure which provides vibration excitation inflight 

SPECIFICATIONS— Definition of vibration quantity versus frequency, usually associated 
with programmatic requirements 


TAP TEST-Measurement of apparent mass or accelerance by tapping on structure with 
small rubber or plastic tipped hammer which incorporates force transducer 

TEST FIXTURE-Adapter hardware which allows test item to be mounted to shaker 

TWO DEGREE-OF-FREEDOM SYSTEM (TDFS) —Vibration model with two masses 


4.0 GENERAL REQUIREMENTS 

4. 1 Criteria for Force Limiting . The purpose of force limiting is to reduce the response of 
the test item at its resonances on the shaker in order to replicate the response at the 
combined system resonances in the flight mounting configuration. Force limiting is most 
useful for structure-like test items which exhibit distinct, lightly damped resonances on the 
shaker. Examples are complete spacecraft, cantilevered structures like telescopes and 
antennas, lightly damped assemblies such as cold stages, fragile optical components, and 
equipment with pronounced fundamental modes such as a rigid structure with flexible feet. 
The amount of relief available from force limiting is greatest when the structural impedance 
(effective mass) of the test item is equal to, or greater than, that of the mounting structure. 
However, it is recommended that notches deeper than 14 dB not be implemented without 
appropriate peer review. Force limiting is most beneficial when the penalties of an artificial 
test failure are high. Sometimes this is after an initial test failure in a screening type of test. 

4.2 Instrumentation . 

4-2.1 Piezo-electric Force Gages . The use of piezo-electric force gages for force limiting is 
highly recommended over other types of force measurement means such as strain gages, 
armature current., etc. The advent of piezo-electric, quartz force gages has made the 
measurement of force in vibration tests almost as convenient and accurate as the 


measurement of acceleration. The high degree of linearity, dynamic range, rigidity, and 
stability of quartz make it an excellent transducer material for both accelerometers and force 
gages. Similar signal processing, charge amplifiers and voltage amplifiers, may be used for 
piezo-electric force gages and accelerometers. However, there are several important 
differences between these two types of measurement. Force gages must be inserted 
between (in series with) the test item and shaker and therefore require special fixturing, 
whereas accelerometers are placed upon ( in parallel with) the test item or shaker. The total 
force into the test item from several gages placed at each shaker attachment may be obtained 
by simply using a junction to add the charges before they are converted to voltage, whereas 
the output of several accelerometers is typically averaged rather than summed. Finally, 
piezo-electric force gages tend to put out more charge than piezo-electric accelerometers 
because the force gage crystals experience higher loading forces, so sometimes it is 
necessary to use a charge attenuator before the charge amplifier. 

4.2.2 Force Gage Preload . Piezo-electric force gages must be preloaded so that the 
transducer always operates in compression. Having a high preload and smooth transducer 
and mating surfaces minimizes several common types of gage measurement errors, e.g. 
bending moments being falsely sensed as tension/compression. However, using flight 
hardware and fasteners, it is usually impossible to achieve the manufacturers recommended 
preload. In addition, sometimes it is necessary to trade-off transducer preload and dynamic 
load, particularly moment, carrying capability. The two requirements for selecting the 
preload are that it is sufficient to prevent unloading due to the dynamic forces and moments 
and that the maximum stress on the transducers does not exceed that associated with the 
manufacturer’s recommended maximum load configuration. 

Transducer preloading is applied using a threaded bolt or stud which passes through the 
inside diameter of the transducer. With this installation, the bolt or stud acts to shunt past 
the transducer a small portion of any subsequently applied load, thereby effectively 
reducing the transducer’s sensitivity. Calibration data for the installed transducers is 
available from the manufacturer it they are installed with the manufacturer’s standard 
mounting hardware. Otherwise, the transducers must be calibrated in situ. 

4.2.3 Force Gage Calibration. The force gage manufacturer provides a nominal calibration 
for each transducer, but the sensitivity of installed units must be determined in situ, as 
discussed in the previous paragraph. This may be accomplished either quasi-statically or 
dynamically. Using the transducer manufacturer’s charge amplifiers and a low noise cable, 
the transducers will hold their charge for several hours, so it is possible to calibrate them 
statically with weights or with a hydraulic loading machine. It is recommended that the 
calibration be performed by loading the transducers, zeroing out the charge, and then 
removing the load, in order to minimize the transient overshoot. 

The simplest method of calibrating the transducers for a force limited vibration test is to 
conduct a preliminary low level sine sweep or random run and to compare the apparent 
mass (ratio of total force in the shaker direction to the input acceleration) measured at 
frequencies much lower than the first resonance frequency with the total mass of the test 
item. Typically the measured force will be approximately 80 to 90 % of the weight in the 
axial direction and 90 to 95% of the weight in the lateral directions, where the preloading 
bolts are in bending rather than in tension or compression. Alternately, the calibration 
correction factor due to the transducer preloading bolt load path may be calculated by 
partitioning the load through the two parallel load paths according to their stiffness; the 
transducer stiffness is provided by the manufacturer, and the preload bolt stiffness in 
tension and compression or bending must be calculated. 
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4.2.4 Force Gage Combinations It is recommended that the total force in the shaker 
excitation direction be measured in a force limited vibration test. The total force from a 
number of gages in parallel is readily obtained using a junction box which sums the 
charges, and therefore the forces, before conditioning the signal with a charge amplifier. 
An alternative is to specify limits for the force at individual attachment positions, but this is 
not recommended. Since vibration tests are normally conducted sequentially in three 
perpendicular axes, it is convenient to employ triaxial force transducers. Sometimes it is 
necessary to limit the cross-axis force and the moments in addition to the in-axis force; this 
is particularly the case in tests of large eccentric test items such as spacecraft. For these 
applications, the six force resultant forces and moments for a single node may be measured 
with a combination, commonly four, of triaxial force transducers and a voltage summer. 

4.2.5 Accelerometers . Accelerometers on the fixture are also required in force limited 
vibration tests in order to control the acceleration input to the acceleration specification at 
frequencies other than at the test item resonances. In addition, it is often convenient to use a 
limited number of accelerometers to measure the response at critical positions on the test 
item. These response accelerometers may be used only for monitoring or, if justified by 
appropriate rationale, for response limiting in addition to the force limiting. 

4.3 Fixturing . The preferred method of configuring the force gages is to sandwich one 
gage between the test item and conventional test fixture at each attachment position and use 
fasteners which are longer than the conventional ones to accommodate the height of the 
gages. In this configuration, there is no fixture weight above the transducers and the gage 
force is identical to the force into the test item. Sometimes the preferred approach is 
impractical, e. g. if there are too many attachment points or the attachments involve shear 
pins in addition to bolts. In these cases it may be necessary to use one or more adapter 
plates to interface the transducers to the test item. The requirement is that the total weight of 
the adapter plates above the force gages does not exceed ten percent of the weight of the test 
item. This limitation is necessary because the force gages read the sum of the force required 
to accelerate the interface plate and that delivered to the test item. If the fixture weight 
exceeds the 10% criterion, force limiting will only be useful for the first one or two modes 
in each axis. Use of a circuit to subtract the interface plate force in real time, is not 
recommended because of the errors that result when the interface plate is not rigid. The use 
of armature current to measure shaker force is also not generally useful, because the weight 
of the armature and fixturing typically are much greater than 10 % of that of the test item. 

4.4 Force Specifications . Force limits are analogous and complementary to the acceleration 
specifications used in conventional vibration testing. Just as the acceleration specification is 
the frequency spectrum envelope of the inflight acceleration at the interface between the test 
item and flight mounting structure, the force limit is the envelope of the inflight force at the 
interface. In force limited vibration tests, both the acceleration and force specifications are 
needed, and the force specification is proportional to the acceleration specification. 
Therefore force limiting does not compensate for errors in the development of the 
acceleration specification , e.g. undue conservatism or lack thereof. These errors will carry 
over into the force specification. Since inflight vibratory force data are lacking, force limits 
are usually derived from coupled system analyses and impedance information obtained 
from measurements or finite element models (FEM). Fortunately, considerable data on the 
interface force between spacecraft and components are becoming available from spacecraft 
acoustic tests, and semi-empirical methods of predicting force limits are being developed. 

4.4.1 Analytical Force Limits . Analytical models and methods of obtaining impedance 
information to use in these models are discussed in Section 5.0 Detailed Requirements. 
Here, the general requirements for analytical force limits are discussed. It is required that 
analytical models used to predict force limits take into account the resonant behavior of both 
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the source (mounting structure) and the load (test item) and that the models incorporate 
impedance information, data or FEM, on both the source and the load. The models 
discussed in the Detailed Requirements section are two-degree-of-freedom system (TDFS) 
models, in which the coupled source and load are each described by a single resonant 
mode. In more complex models, the source and load may have many modes. In the early 
stages of a program, before hardware exists, strength of materials or FEM models are often 
used to determine the modal parameters of the source and load. Later in the program, 
before the vibration tests of flight hardware, it is recommended that the modal parameters 
be updated with impedance data measured in tap tests on the mounting structure and in the 
shaker tests of the test item. The coupled source and load models are exercised with some 
representative excitation of the source, and the envelope (or peak values) of the interface 
acceleration and interface force frequency response functions (FRF) are calculated, 
preferably in one-third octave bands. Finally, the ratio of the interface force envelope to the 
acceleration envelope is calculated from the model and the force limit specification is 
calculated by multiplying the conventional acceleration specification by this ratio. (It is 
essential that the ratio of envelopes or peaks, not of the FRF’s, be calculated.) 

4.4.2 Semi-empirical Force Limits . The alternative semi-empirical approach to deriving 
force limits is based on the extrapolation of interface force data for similar mounting 
structure and test items. A general form for a semi-empirical force limit for sine or transient 
tests is from “Taming the General-Purpose Vibration Test”: 


F, = CM„A, (la) 

where F, is the amplitude of the force limit, C is a frequency dependent constant which 
depends on the configuration, M 0 is the total mass of the load (test item), and A s is the 
amplitude of the acceleration specification. The form of Eq. la appropriate for random 
vibration tests is: 


S ff =C 2 M 0 2 S aa (lb) 

where S FF is the force spectral density and S AA the acceleration spectral density. 

As shown in “Verification of Force and Acceleration Specifications for Random Vibration 
Tests of Cassini Spacecraft Equipment”, interface force data measured between three 
instruments, each weighing approximately 60 lb., and JPL’s Cassini spacecraft in acoustic 
tests of the development test model spacecraft fit Eq. lb with C equal to unity at 
frequencies up to and including the fundamental resonance of the test item and then with C 
rolling off as one over frequency at higher frequencies. It is required to show similarity 
between the subject hardware configuration and the reference data case or to justify the 
scaling used for any extrapolation, in order to use semi-empirical force limits. 

4.4.3 Quasi-static Design Verification . The quasi-static design of aerospace components is 
often based on a specified acceleration of the center-of-gravity (C.G.) of the component. 
However, the C.G. of a flexible body is a virtual (not a real) point and its acceleration 
cannot be accurately measured with an accelerometer in a vibration test, particularly at 
frequencies above the fundamental resonance. However, Eq. la with C equal to unity and 
A s equal to the C.G. acceleration is Newton’s 2nd law. Thus limiting the external force to 
the product of total mass times the quasi-static design limit, or some fraction thereof, is the 
recommended method of validating quasi-static designs in vibration tests. 
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4.5 Control System . Most of the current generation of vibration test controllers have the 
two capabilities needed to implement force limiting. First, the controller must be capable of 
extremal control, sometimes called maximum or peak control by different vendors. In 
extremal control, the largest of a set of signals is limited to the reference spectrum. (This is 
in contrast to the average control mode in which the average of a set of signals is compared 
to the reference signal.) Most controllers used in aerospace testing laboratories support the 
extremal control mode. The second capability required is that the controller must support 
different reference spectra for the response limiting channels, so that the force signals may 
have limit criteria specified as a function of frequency. Controllers which support different 
reference spectra for limit channels are now available from most venders and in addition 
upgrade packages are available to retrofit some of the older controllers for this capability. If 
the controller does not have these capabilities, notching of the acceleration specification to 
limit the measured force to the force specification must be done manually in low level runs. 


4.6 Test Planning Considerations . Several considerations need to be addressed in the test 
planning if force limiting is to be employed. First the size, number, and availability of the 
force transducers need to be identified as well as any special fixturing requirements to 
accommodate the transducers. Next, the approach for deriving and updating the force 
specification needs to be decided. Finally the control strategy, which in special cases may 
include cross-axis force, moment, individual force, and response limiting in addition to or 
in lieu of the in-axis force, must be decided and written into the test plan. In some 
instances, the control strategy will be limited by the control system capabilities. In all cases, 
it is recommended that the control strategy be kept as simple as possible, in order to 
expedite the test and to minimize the possibility of mistakes. 


5.0 DETAILED REQUIREMENTS 

5.1 Derivation of Force Limits . As the force limiting technology matures, there will 
eventually be as many methods of deriving force limits as there are of deriving acceleration 
specifications. Herein several acceptable methods are described. 

Force spectra have typically been developed in one-third octave bands (see example in 
Section 6.2), but other bandwidths, e.g. octave or one-tenth octave bands, may also be 
used. Force limiting may usually be restricted to an upper frequency encompassing 
approximately the first three modes in each axis; which might be approximately 100 Hz for 
a large spacecraft, 500 Hz for an instrument, or 2000 Hz for a small component. It is 
important to take into account that the test item resonances on the shaker occur at 
considerably higher frequencies than in flight. Therefore care must be taken not to roll off 
the force specification at a frequency lower than the fundamental resonance on the shaker 
and not to roll off the specification too steeply, i.e. it is recommended that the roll-offs of 
the force spectrum be limited to approximately 9 dB/octave. 

5.1.1 Simple TDFS . The simple Two-Degree-of-Freedom System (TDFS) method of 
deriving force limits is described in “Vibration-Test Force Limits Derived from Frequency- 
Shift Method”. The basic model is shown in Figure 1 . The model represents one vibration 
mode of the source (system 1) coupled with one vibration mode of the load (system 2). 
Figure 2 shows the ratio of the interface force spectral density S FF to the input acceleration 
spectral density S AA , normalized by the load mass M 2 squared, as a function of the mass 
ratio M/M,, calculated from the simple TDFS. When this mass ratio is very small, there is 
no force limiting effect; the force spectral density asymptote is the load mass M 2 squared 
times the input acceleration spectral density times the quality factor Q, squared. The ratio of 
this asymptotic value of the force to the force limit at larger values of M/M,, is the 
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6.2 Force Specification Example . Appendix C is a spread sheet calculation of the force 
specification for an instrument (CRIS) mounted on a honeycomb panel ot a spacecraft 
(ACE) using three methods: the simple TDFS, the complex TDFS, and the semi-empirical. 


6.3 Definition of Symbols 

A = interface acceleration 

Ab = base acceleration 

Ao = free acceleration of source 

As = acceleration specification 

C = dashpot constant 

C = constant 

F = interface force 

Fs = force specification or limit 

k = spring stiffness 

k = physical stiffness matrix 

M 0 = total mass 

M = residual mass 

m = modal mass 

M = apparent mass, F/A 

m = physical mass matrix 

M = modal mass matrix 

Q = dynamic amplification factor 

S AA = acceleration spectral density 

S FF = force spectral density 

u = absolute displacements 

U = generalized modal displacement 

(J) = mode shape 

to = radian frequency 

to o = natural frequency of uncoupled oscillator 

Subscripts 

1 — source oscillator 

2 = load oscillator 

F = unrestrained (free) 

P = prescribed 

N = modal set 

R = rigid body set 

n = single mode 

p = reaction force direction 

q = prescribed acceleration direction 
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FIGURE 1 . Simple TDFS of Coupled Oscillators 
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FIGURE 2. Normalized Force Spectrum for Simple TDFS 
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0.25. 0 5 

103 

99 

86 

63 

38 

22 

14 

12 

! 1 

0.25. 0.25 

28 

28 

27 

23 

16 

10 

8 

7 

7 

0.25. 0.125 

8 

8 

8 

8 

7 

5 

5 

4 

4 

TABLES I. Normalized Force Spectrum for Complex TDFS with 0=20 
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APPENDICES 

Appendix A-Equations for Calculating the Simple TDFS Force Limits 

The force limit is calculated for the TDFS in Fig. 1 with different masses of the source 
and the load oscillators. For this TDFS, the maximum response of the load and therefore 
the maximum interface force occur when the uncoupled resonance frequency ot the load 
equals that of the source. For this case, the characteristic equation is that of a classical 
dynamic absorber, from “Vibration-Test Force Limits Derived from Frequency-Shift 
Method”: 

(co/co o ) 2 = 1+ (m 2 /m,)/2 ± [(m/m,) + (m 2 /m,)"/4)] 0 (Al) 

where to, is the natural frequency of one of the uncoupled oscillators, m, is the mass of the 
source oscillator, and m 2 is the mass of the load oscillator in Fig. 1. The ratio of the 
interface force S KF to acceleration S AA spectral densities, divided by the magnitude square 
of the load dynamic mass m 2 , is: 

S hT /(S AA m 2 2 ) = [1+ (to/w 0 ) 2 /Q 2 2 ] /{ [1- (o)/to 0 ) 2 ] 2 + (0 )/g> o ) 2 /Q 2 2 } (A2) 

where is the quality factor, one over twice the critical damping ratio, of the load. 

The force spectral density, normalized by the load mass squared and by the acceleration 
spectral density, at the two coupled system resonances is obtained by combining Eqs. (Al) 
and (A2) For this TDFS the normalized force is just slightly larger at the lower resonance 
frequency of Eq. (Al). The maximum normalized force spectral density, obtained by 
evaluating Eq. (A2) at the lower resonance frequency, is plotted against the ratio of load to 
source mass for three values of Q 2 in Fig. 2. 


Appendix B--Calculation of Effective Mass 


Applying the rationale of “Equivalent Spring-Mass System: A Physical Interpretation and 
subdividing the displacement vector into unrestrained absolute displacements u, and 
prescribed absolute displacements u p , the equilibrium equation is: 


Let: 


[ nip F 

1 nip P ] 

d 2 Up. /dt 2 
f \ + 

[ kj; F 1 k FP ] 

[ m P p. 

l FHpp 1 

d 2 u p /dt 2 

[ kpj; 1 kpp ] 



[ 4>n 1 4> P ] 

U N 

{u} = 

<)> U = 


{-} 



[0 1 I PP ] 

Up 


u F 

{-} 

u P 


fr 

{-> 

fr 


(Bl) 


(B2) 


Where <|) N are normal modes and (|) R are rigid body modes associated with a kinematic set of 
unit prescribed motions, and U N is the generalized modal relative displacement and U p is 

the generalized prescribed absolute displacement. Substituting and pre-multiplying by 0 T 
yields: 


15 



where: 


[M nn 1 

m np ] 

d 2 u N /dr 

[ m t np 

1 Mpp ] 

1 

d 2 Up/dt 2 


[ co N M nn 
[ 0 


0 ] 
0 ] 



F k 

{— -} (B3) 

Fp 


^NN — 0N m FF 0N (B4) 

M np = <(> n t m FF <t> p + 0 n t m, p I pp (B5) 

M pp — I pp m PI> I pp + I p p rripp 4> F + <J> P m, ;p I pp + (J> p r m F(; <j) p (B6) 

F P = Ipp fp (B7) 

For: d 2 U F /dr = U p =F F = 0,d 2 U fl /dt 2 = -6) 2 U n ,andforU„ = I: 

M n p ] = - Fp / co n 2 (B8) 


where n indicates a single mode. (Note that M nP r is in mass units.) M nP is sometimes called 
the elastic-rigid coupling or the modal participation factor for the nth mode. If the 
model is restrained at a single point, the reaction (F ) in (B8) is the SPCFORCE at that 
point in a NASTRAN modal analysis. 

The initial value ot M pp is the rigid body mass matrix. If a Gaussian decomposition of the 
total modal mass in (B3) is performed, it subtracts the contribution of each normal mode, 
called the effective mass: 

M np‘ M nn~' M nP ’ (B9) 

from Mpp", which is the residual mass after excluding the mass associated with the 
already processed n modes. 

Consider the ratio of the reaction force in a particular direction p, to the prescribed 
acceleration in a particular direction q; the effective mass, M n T M nn 1 M n , is the same as the 
contribution of the n th mode to this ratio, divided by the single-degree-of-freedom 
frequency response factor. The sum of the common-direction effective masses for all 
modes is equal to the total mass, or moment of inertia for that direction. Values of the 
effective mass are independent of the modal normalization. 
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Appendix C— Force Specification Example 
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